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Abstract 


This thesis presents a state space approach to optimal estimation of channel processes, such 
as inphase, quadrature, square envelope, phase, etc., of wireless fading channels. The mod¬ 
els are derived from the Doppler power spectral density using factorization and realization 
techniques. Several simulations are performed when the channel is Rayleigh and Ricean dis¬ 
tributed under both flat and frequency-selective fading assumptions. The simulation studies 
illustrate the tracking properties of the inphase, quadrature and square envelope estimators 
for various received signal-to-noise ratios. 
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Chapter 1 


Introduction 


1.1 Mobile Radio Channel 


Three fundamental components of a wireless radio communication system are the trans¬ 
mitter, channel and the receiver. The transmitter sends messages to the receiver through 
a channel. These messages, called the transmitted signals, are some kinds of electromag¬ 
netic waves, which obey the well known Maxwell’s equations while propagating through the 
channel. The physical channel is usually the atmosphere (the free space) in wireless radio 
communications. One of the most important characteristics in wireless radio communications 
is fading, which causes the received signal power to change with time. The change of the 
power at some place or on some time interval is so big that the transmitted signal vanishes 
into the void. A concrete example is signal propagation in an urban area, where natural 
and man-made objects are concentrated. This kind of signal propagation gives rise to the so 
called multipath fading, in which the received electromagnetic waves consist of superposition 
of a number of the transmitted signals, which each one has a differential amplitude and time 
delay. 

Statistical techniques have been introduced to describe this kind of channels. In view of 
this, the channel statistical properties, like mean, variance, autocorrelation, power spectral 
density and so on, have been investigated. 

Three basic mechanisms which affect signal propagation in wireless radio communications 
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are reflection 1 , scattering 2 and diffraction 3 . When the electromagnetic waves propagate 
through the atmosphere, they undergo reflections, scattering from various surfaces of objects 
and diffractions. These phenomena can be described by stochastic variables for developing 
various mathematical fading channel models. After a channel model is constructed, how 
’’good” the model is has to be we evaluated against real data. It is convenient to simulate 
the channel model in a software simulator at first in stead of implementing the channel in 
hardware. Doing this way, we will save both time and money. 

Channel estimations allow the receiver to approximate the impulse response of the channel 
and explain the behavior of the channel. The knowledge of the channel’s behavior is well- 
utilized in modern radio communications. It must be stressed here that channel estimation is 
only a mathematical representation of what is truly happening. A ’’good” channel estimation 
is one where some sort of error minimization criteria is satisfied. Kalman filtering theory 
which belongs to the least-square optimal estimation techniques is used in this paper to 
estimate the constructed channel inphase, quadrature and square envelope. 

Emerging in the area of systems and control theory, the Kalman filter provides a method 
for generating an optimal estimation of the system state. That is, given a signal model 
that consists of a linear dynamical system driven by stochastic white noise processes, the 
Kalman filter provides an optimal way of extracting a signal from noise by exploiting a state 
space signal model. The prerequisite of using Kalman filter is the existence of a state space 
channel model. The approximation, factorization and realization techniques will help us for 
constructing the state space channel model. 

In general, the multipath propagation of the transmitted signal can be divided into large- 
scale fading 4 and small-scale fading 5 . Two typical small-scale fading channels are Rayleigh 
and Rican, which we will introduce in the following section. 

1 Reflection occurs when a propagating electromagnetic wave strikes a smooth surface with very large 
dimensions compared to the RF signal wavelength [31]. 

Scattering occurs when a radio wave impinges on either a large rough surface or any surface whose 
dimensions are on the order of the signal’s wavelength or less, causing the reflected energy to spread out 
(scatter) in all directions [31]. 

3 Diffraction occurs when the radio path between the transmitter and receiver is obstructed by a dense 
body with large dimensions compared to the wavelength, causing secondary waves to be formed behind the 
obstructing body. Diffraction is a phenomenon that accounts for RF energy traveling from transmitter to 
receiver without using a line of site path between these two [31]. 

4 Large-scale fading: represents an average received signal power attenuation or path loss due to motion 
over large areas [31]. 

5 Small-scale fading: refers to the dramatic changes in received signal amplitude and phase that can be 
experienced as a result of small changes in spatial separation between a transmitter and receiver [31]. 
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1.2 Aulin’s Channel Model 

(i) Representation. In small-scale fading statistics, several multipath fading channel 
models were reported in the literatures ([1], [2], [3]). Ossanna’s model [1] is one of first, and 
assumes existence of a direct path between the transmitter and receiver. Clarke’s model [2], 
which includes Ossanna’s model as a special case, has the property of signal scattered in the 
vicinity of the receiver. Aulin’s model [3], is a generalization of the Clarke’s model, because 
it considers a three dimensional wave propagation model. Our investigation in this thesis, is 
based on the Aulin three dimensional model which is shown in Figure 1.2.1. 


z 



Figure 1.2.1: A wave component in three dimensions 

It is important to note that additional models can also be included, such as those in Gans 
[5]. It assumed that at any point of the received field, the total wave consists of superposition 
of'lV plane waves, which each one travels via a different path. Let 0(0,0,0) denote the zero 
reference phase point; let the n th incoming plane wave, traveling through the n th path, be 
denoted by E n (t) with its parameters {a n , j3 ni (f> n , c n }^ =1 . Here c n and 4> n are the amplitude 
and phase of the n th plane wave. Suppose the transmitted signal is Re{e^ ct }, then at point 
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O, the n th incoming wave can be represented by 


E n (t) = Re{c n e j(t>n e ju)ct }, 1 < n < N, 

where Re{x } denotes the real part of the complex-valued quantity x. Further, at any receiv¬ 
ing point O'(xo,yo, z 0 ), the n th plane wave is 6 

E n (t) = Re {cne^e?^ 1 -^} 

r 27T 

= C n cos \ijj c t - — (xq cos a n cos (3 n + y 0 sin a n cos f3 n + z 0 sin /3 n ) + <p n , 

where A is the wavelength. Next, let the point O ! (x 0 , y 0 , z 0 ) move in the received field with 
a velocity v on the x — y plane, which has a direction with an angle 7 to the x — z plane, 
then 0'(x 0 , y 0 , z 0 ) moves to a new position 0"(x 0 + vt cos 7 , y 0 + vt sin 7 , zq) at time t. Thus, 


E n ( t ) = c n cos (u n t + 0 n ) cos 0 o c t - c n sin (u n t + 0 n ) sin uj c t 
— I n (t ) cos u c t - Q n (t) sin a; c t, 

= Re{c n e>* n ei“° t }, 


where 


OJn 


2-kv 

~ 


cos (7 - a n ) cos p n , 0 n - 


2wzo 


sin P n 


2n^xl + 


Vo 


X 


cos(a n + £) cos p n + 4> n . 


Here u n is the Doppler shift, 0 n is the phase associated with the n th wave, I n (t), Q n (t) are the 
inphase and quadrature components, c n = \Jln(t) + Q^(f), $n(t) = tan ~ l {Q n {t)/I n {t)) — 
oj n t + On and £ = tan _ 1 (y 0 /^o)- Thus, at any receiving point, the superposition of the 
incoming waves is given by 

m = EKm 

n =1 

N N 

= £ Cn cos (u n t + On) cos Uj c t - ^2 Cn sin (u n t + 0 n ) sin u c t 

n=l n=1 

= I(t) cos u) c t — Q(t ) sinu; c f, 


where /(f) = En=i c n cos (u n t+O n ), Q(t) = En=i °n sm(uj n t+Q n ). When N is sufficiently large 
(typically N > 6 ), both /(t) and Q(t) have the statistical property of Gaussian distribution, 
and so does E(t). 


6 see Appendix A. 1.1 for details 
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(ii) Mean, Variance and Autocorrelation. It is easily verified that 7 

E[I(t)} = E[Q(t)} = E[E(t)} = 0, 

E[E(t)E(t + r)] = a(r ) coso; c r — c(r) sinu> c T, 


where 


a(r) = £[/(£)/(£ + r)] = E\Q(t)Q(t + r)] = YULi -E[coso; n r], 
c(r) = £?[/(i)Q(i + r)] = E[l(t + r)Q(t)} = § Eli £[smw n r]. 


Here Eq is a positive constant. From above equations, we can easily get the second moment 
of I(t ) or Q(t) by letting r = 0, that is, a(0) = Eq/2. Further, assuming that a n is uniformly 
distributed on [0,27r], then we have 

o(r) = y j J a ( 2 -^ cos 0)pa ff))dp, (1.2.1) 

c(t) = 0, (1.2.2) 


where J 0 (-) is a Bessel function of the first kind of zero order. Thus, we conclude that I(t) 
and Q(t ) are iid with density iV(0; Eq/2). 

If there is a direct path between the transmitter and receiver, the means of the inphase 
and quadrature are finite and not zero, especially, define 


where 


H — E[I d (t)\ = r 0 cos(ta 0 £ + 9 0 ), v — E[Q d (t)] = r 0 sin(ta 0 £ + 6 0 ), 


2iw . . ^ 2ttz 0 . „ . 

w 0 = —— cos (7 - a 0 ) cos /?o, #0 =- 7 — sm/3 0 + 4>o 


with {r o ,w o , 0 o} are unknown constants, then 


E(t) = (/(£) + I d (t )) cost o c t - ( Q(t ) + Q d (t )) sinw c £. 

Obviously, the inphase and quadrature components are still Gaussian with densities iV(/j; o(0)) 
and N(u-,a( 0)). Usually, the direct path between the transmitter and receiver is called the 
specular component or the line of sight (LOS). 

(iii) Received Signal Amplitude Distributions. The received signal amplitude 
which is represented by the inphase and quadrature components has the following distribu¬ 
tions. 


7 see Appendix A.1.2 for details 
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Rayleigh Distribution. Since I(t), Q(t ) ~ lV(0;a(0)). It is easily to prove that the signal 
envelope r = \JI 2 {t) + Q 2 (t) is Rayleigh distributed with pdf 8 


f( r ) = ■ r . e -r 2 /2a(0) r > Q 
J[) o( 0) ’ 

and £[r] = y®, Var(r) - (2 - f)o(0). 

Ricean Distribution. Since /(f) + /d(f) ~ N(fj,\a( 0)), <5(t) + (3</£) ~ lV(i/;a(0)). Obvi¬ 
ously, the signal envelope r = ^(/(f) + I d (t )) 2 + (Q(t) + Qd(t)) 2 is obeyed Ricean distribu¬ 
tion with pdf 9 


f( r ) = 


°( 0 ) 


2,2 

-. r - + r a r 
e 2 ° (°) Jo 



r > 0, 


where / 0 (-) is the modified Bessel function of the first kind of zero order, and E[r 2 ] = 
Tq + 2a(0), Var(r 2 ) = 4a 2 (0) — 4a(0)rQ. 

(iv) Power Spectrum. Suppose that pp((3) = 5 (ft) which means the incoming waves 
travel only on x — y plane. Integrate (1.2.1), we get ao(r) = ^ Joi^j 1 - cos/?). Thus, the 
power spectral density function of lit ) or Q(^) is obtained by taking the Fourier transform 
to a 0 (r), 

$>(/) = ^{ao(T)} 


' £n l A 1 

< 2 27r VM£) 2 ’ 


I/I < !; 

otherwise. 


The processes are strictly band-limited within the maximum Doppler frequency shift fa (see 
Figure 1.2.2). 

Suppose that 


101 <IW<|; 

0, otherwise, 

which means the incoming waves travel in the three dimensional space. Take the Fourier 
transform to (1.2.1), then we obtain the Doppler power spectral density function (DPSD) 
for small angle /? m 10 , 



8 see Appendix A. 1.3 for details 
9 see Appendix A. 1.4 for details 
10 see Appendix A. 1.2 for details 
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f = 910 MHz, v= 40 km/h, f. = 33.704 Hz 
c d 



Figure 1.2.2: The Doppler power spectrum of pp(/3) — S{P) 


Sn(f) = FMt)} 


t 


0, 



Eq 1 A 
2 2 sin v ’ 


j COS 1J,„ 

< I/I < x; 


E a 1 A 
2 27rsin/3 m v 


— arcsin 


2cos 2 /? m —1—(^7) 

i-(^r ) 2 



I/I < 5cos/? m 


Let the carrier frequency / c = 910 MHz, v change from 5 km/h to 200 km/h, then we 
obtain a series of the Doppler power spectrums shown on Figure 1.2.3. 


1.3 Concepts about Fading Channel 

In this section, we are going to introduce some concepts about fading channel for the 
purposes of simulation of the channel estimation. Let the transmitter and/or the receiver 
move, then there is a phase shift, called Doppler shift, on the received signal. The value 
of Doppler shift depends on the velocity of the moving object, the carrier frequency of the 
transmitted signal and the angle between the direction of the object’s movement and the 
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Frequency Hz 


Figure 1.2.3: 3-D power spectrums 

direction of the signal arrival. The maximum Doppler shift is defined by fd = v/X, which is 
a function of the velocity and the length associate with the carrier wave. 

Bd, the Doppler spread, is a measure of the spectral broadening caused by the time rate 
of change of the mobile radio channel and is defined as the range of frequencies over which 
the received Doppler spectrum is essentially non-zero. That is, Bd = 2 fd- 

The reciprocal of the Doppler spread is defined as the channel coherence time, T c = 
1 /Bd, which is the temporal interval over which the channel characteristics remain relatively 
constant. The time-variant nature of the channel can be viewed in terms of fast fading and 
slow fading. 

A channel is called fast fading when T s » T c , where T s is the time duration of a 
transmission symbol (signaling interval). It describes a condition where the time duration 
in which the channel behaves in correlated manner is short compared to the time duration 
of the symbol. Therefore, it can be expected that the fading character of the channel will 
change several times while a symbol is propagating. A channel is referred to slow fading when 
T a « T c . That is, the time duration which the channel behaves in a correlated manner 
is long compared to the time duration of a transmission symbol. Thus, one can expect the 
channel state to remain unchanged during the time in which a symbol is transmitted. 
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B c , the channel coherence bandwidth, is a statistical measure of the range of frequencies 
over which the channel passes all spectral components with approximately equal gain and 
linear phase. That is, a signal’s spectral components in that range are affected by the channel 
in a similar manner. It is defined as B c = 1/T m . 

T m , called the multipath spread or the delay spread of the channel, is the standard 
deviation (or root-mean-square) value of the delay of reflections, weighted proportional to 
the energy in the reflected waves. In a fading channel, the relationship between the delay 
spread T m and the symbol duration T s can be viewed in terms of two different categories, 
frequency-selective fading and frequency-nonselective or flat fading. 

A channel is said to exhibit frequency-selective fading if T s <T m , or when the bandwidth 
of the transmitted signal is less than the bandwidth of the channel, B s < B c . A channel 
is said to exhibit frequency-nonselective fading if T s > T m , or when the bandwidth of the 
transmitted signal is greater than the bandwidth of the channel, B s > B c . 

1.4 Outline 

A key to our approach is the Doppler power spectral density function which we introduced 
in section 1.2. Based on it, we will do the following works in this thesis: 

1. Approximation of the Doppler power spectral density function for obtaining a system 
transfer function. 

2. Factorization of the system transfer function as in [11]. 

3. Realization of the system transfer function for getting a state space system model. 

4. Estimation of the state space system model. 

5. Simulation of the state space system estimations. 

The thesis is presented as follows. 

" Chapter 2, introduces the multipath fading channel model, and some preliminary material 
from random signal and systems. Next, a state space model is derive from Doppler power 
spectral density, using approximation, factorization and realization techniques. Then, we 
simulate some characteristics of Rayleigh and Ricean state space fading channel models. 
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Chapter 3, introduces the state space channel estimation problem. It presents the Kalman 
filtering algorithms for estimating the inphase, quadrature and square envelope. In order to 
use computer to perform these estimates, we discretize the model. In addition, conditional 
distributions for the channel envelope and phase are derived in closed form. 

Chapter 4, presents several case studies of the estimation problem considered in chapter 
3, assuming that the state space channel model parameters {A,B,G,D} are known. Simu¬ 
lations include inphase, quadrature and square envelope estimates for Rayleigh and Ricean 
channels. Both fiat slow fading and flat fast fading are considered. Finally, a multipath 
state space fading channel model which involves one Ricean and five Rayleigh’s components 
is simulated. 

Chapter 5, presents certain topics of future research which emerge from the finding of 
the thesis. 



Chapter 2 

Channel Modeling 


In general, in wireless radio communications, radio propagation is mainly by way of 
scattering from the surfaces of the objects and by diffraction over and/or around them. 
These natural phenomena produce very complicated signal, which is a summation of the 
transmitted signal, at the receiver. Each term of the summation stands for a multipath 
component having an attenuation and a time delay. The number of multipath components, 
attenuation and time delay are all stochastic processes. Based on the descriptions in the 
previous chapter, we will do the following tasks in this chapter: 


1. Representation of a multipath fading channel model. First giving the received signal 
model in lowpass and bandpass form, then we introduce the related background about 
a wide-sense stationary signal which is transmitted through a linear time-invariant 
system. 

2. Factorization of the Doppler power spectral density-(DPSD). We approximate the 
DPSD by a fourth-order rational function which is similar to [6]. 

3. Realization of DPSD by state space models. We present a fourth-order approximation 
of the DPSD and its equivalent state space realization as in [8]. 

4. Simulation of the state space channel models. The inphase, quadrature, envelope 
and phase of the wireless state space Rayleigh and Ricean fading channel models are 
simulated by using Simulink. 


11 
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2.1 Multipath Fading Channel Model 

Suppose we transmit a lowpass signal {si(t )} t >o with a carrier frequency f c through a 
wireless multipath fading channel. The received signal in the lowpass representation is 

N(t) 

Vi(t ) = £ - n ) + 

i=l 

where (tu* (£)}<>o is a complex-valued white Gaussian noise process, iV(£) is the total num¬ 
ber of multipath components at time t, Ti(t,Ti ) is the real amplitude of the I th multipath 
component at time t, <f>j(£, Tj) is the phase shift due to free space propagation of the i th 
multipath component plus any additional phase shifts which are encountered in the channel 
at time t and Tj is the excess time delay on the i th path. Here Tj(£, is called the 

attenuation factor for the received signal on the i th path, which can be used to describe the 
wireless multipath fading channel. Another representation is obtained by using the inphase 
and quadrature components of the i th path. That is, 




= h(t,Ti) +jQi{t,Ti), 


where Ii(t,r) and Qi(t,r) are the inphase and quadrature components of the channel corre¬ 
sponding to the i th path. Equivalently, the envelope representation is 






]/ii(t,n) + Qi(t,Ti)e- 7 tan 



Alternatively, the received signal in the bandpass representation is 

y(t) = i?e| [ J2 r ’i(£» r i) eJ$ ‘ (t,T<)s i( t - R ) +Wi(tj\e jwct 

Consider the lowpass equivalent representation of the transmitted signal Si(t ) = Si(t)+jSQ(t) 
with the inphase and quadrature components Si(t) and SQ(t). Then, the bandpass received 
signal in the envelope format is given by 


y(t) = 


IS(t) f 

iz | y/lifoR) + QKt>Ti)y/Sf(t - nj + S%{t - n ) 

X COS (w c t + (£, Tj) + 6(t, Tj) + a(t, Tj)) | 

+ 1 Jv}(t)+V%(t) COS (w c t + T}(tj) , 



2.2. RANDOM SIGNALS AND LINEAR SYSTEMS 
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where 


8 (t, Tj) = tan 1 


Ii{t,Ti ) ’ 


r)(t) = tan 1 


VQ(t) 

Vj(t) ’ 


a(t, Ti) = tan 1 


S Q (t - Tj) 
Sl(t - Ti) ’ 


The bandpass received signal in the quadratures format is given by 


m , 

y(t) = | [■£(*»- r i) - Qi(t> n)S Q {t - Ti)] cos (u c t) 

- Ii(t, Ti)S Q (t - n) + Qi(t, Ti)Sj(t - T^] sin(w c t)| 

+ V!(t) cos (u c t) - v Q (t) sin (u> c t), 


where {vi(t )} t >o and {vQ(t )} t > 0 are two iid white Gaussian noises with density N(0] cr^). 
The simplest bandpass representation of a received signal is 

m 

y(t) = r i(t,Ti) COS 

i —1 

where {s(t)} t >o is a real bandpass transmitted signal, {v(t)} t >o is a bandpass white Gaussian 
noise. The quadratures representation of (2.1.1) is 

y(t) = (h{t, Ti) cos(u c t) - Qi(t, ^) sin(o > c t))s(t - n) 

( 2 . 1 . 2 ) 


(u c t + <hi(t,Ti))s(t - Ti) + v(t), (2.1.1) 


+ vi(t) cos (u c t) - VQ(t) sin(o > c t) 


2.2 Random Signals and Linear Systems 


Consider a wide-sense stationary signal x(t) which passes through a linear time-invariant 
system having impulse response h(t). Suppose h € Li 1 , then the Fourier transform 

/ +oo 

h(t)e~ j2nft dt 

-oo 

is an absolutely convergent integral for each / E (— 00 , 00 ). Define H E Co 2 , Then the 
integral 

is the space of all Lebesgue measurable functions satisfying 

/ +00 

\f(t)\ p dt < 00 , 1 < p < 00 

-OO 

2 Cq is the space of all bounded-continuous functions satisfying 

f(t) —^ 0 as |t| —y 0 
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/ +oo 

H(f)e^ ft df 

-OO 

is called the inverse Fourier transform of H(f) and is not absolutely convergent unless //(/) 
is also in L\. The linear system theory tells us that the convolution product between two 
time function x and h, called the output of the system, is again a wide-sense stationary 
random process, and the definition is given by 

/ +oo 

h(t — r)x{r)dr. 

-OO 

Thus, y{t) = h(t) * x(t) yields 

Y(f) = H(f)X(f), 

for the Fourier transforms of the input X(f), output Y(f) and impulse response H(f). 
Equivalently, in the Laplace transform domain, the following equation, 

Y(s) = H(s)X(s), S=jU , 

hold 3 for the Laplace transforms of the input X(s), output Y (s) and impulse response H(s). 
Further, define the autocorrelations of the input and output signal by 

Rxx(r) = E[x(t + r)x*(t)}, Ry y {r) = E[y(t + r)y*(t)\. 

Obviously, by using linear system theory, the autocorrelation function of y(t) is given by 4 

Ryy(r) = Rxx(r) * (Mt) * h*(-r)). (2.2.3) 

Moreover, define the power spectral densities of the input and output signals by 

S x {f)= / + °° R xx (r)e-^ T dr, S y (f) = [ + °° i^ y (r) e -^dr, -oo < / < oo. 

J —OO J —OO 

3 The relationship between the Fourier and Laplace transform is 

T{x(t)} = C{x{t)}\ 3=ju 


4 see Appendix A.2 for details 
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Prom (2.2.3), the following relationship can be found for the power spectral density functions 
of the input and output signal 5 

Sy(f) =| H(f) I 2 S x (f). (2.2.4) 

If the power spectral density of the input signal satisfies S x (f) = 1, then, from (2.2.4), 
the power spectral density of the output signal equals the power spectrum | H(f) j 2 . For 
a white noise signal which is a Gaussian process, say w(t), with zero-mean variance one, 
one can easily get Syj(f) — l 6 . Compared with the discussions in the previous chapter, 
| H(f ) | 2 is actually the DPSD, and S y (f) is the PSD of the inphase (or quadrature). Thus, 
one can produce the inphase and quadrature branches based on the equation (2.2.4). In 
order to generate the inphase and quadrature components, mathematically, one should find 
an impulse response h(t) which is related to Sr>(f). In the next section, we will introduce 
a rational transfer function H(s) whose square magnitude will approximate to the power 
spectrum density. In other words, we should find a rational transfer function H(s), which 
satisfies 


S D {f) =| H{s) | 2 , s = joj. 


2.3 Factorization of Doppler Power Spectral Density 


Consider factorization of the Doppler power spectral density function 5 d(/). If the Paley- 
Wiener condition L< oo is satisfied, which implies that So(f) is factorizable, 
then we have |#(s)| 2 = So(f), s = ju. Unfortunately, here, the Paley-Wiener condition 
is not satisfied because So(f) is band limited. But if Sr>(f) is approximated by an even 
rational transfer function which is factorizable, and the key properties of the DPSD are not 
affected in this approximation, then this approach will satisfy the Paley-Wiener condition. 
That is, letting 


S D (f) = A 2 


n?=i(s + Qj)(— s + a,i) 
n?=i( s + bi)(—s + bl) 


s=jw 


defining 


H(s ) = B 


nr=i ( 5 + qQ 

W=i(s + bi) 


5 see Appendix A.2 for details 
6 see Appendix A.7.1 for details 
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then 


S D (f) = H(s)H(-s) =| H(s) I 2 , 

where Re(ai) > 0, Re(bi) > 0 and A 2 = \B\ 2 . We choose a second-order rational transfer 
function 

k 

s 2 + 2(u n s + u\ 

to approximate the DPSD. Here u n is the undamped natural frequency, ( is the damping 
ratio and k is the gain. That is, So(f ) is approximated by a fourth-order rational polynomial 
in / through 

S D (f) « H(ju)H(—ju). 


One approximation is to interpolate the points <Sx>(0) and SD(f max ) by determining the 
parameters {u n , (, A;}. In this case, simple algebra 7 shows that 


C = 


1 

1 - A 

^1 

1 


Sp( 0) 

^D^fmax) 


U n = 


2 nf„ 




k^u 2 n y/So{ 0 ), 


(2.3.5) 


where 


fmax fd COs(/3 m ), S[)(fmax) 


Eq 


^f d sin P„ 


(2.3.6) 


Sd(0) = . 7 ~r - arcsin(2cos 2 (3 m - 1)]. (2.3.7) 

47T Jd Sill f Jrn " 

Figure 2.3.1 shows Sp(f ) and its approximation using the above approach for E 0 = 2, 
f c = 910 MHz, Pm — 10°, v — 5 km/h and v = 120 km/h. Figure 2.3.2 shows S D (f ) 
and its approximation using the same approach for E 0 = 2, f c = 910 MHz, P m = 40°, 
v = 5 km/h and v = 120 km/h. Clearly, one can increase the Doppler spread by increasing 
the velocity. On the other hand, one can also reduce the Doppler spread by increasing P m 
while the velocity is fixed. For a much better approximation, one may increase the degree of 
the numerator and denominator of H(s). However, for high order approximations it will be 
difficult to derive explicit equations which are related to the system parameters {u n , </, k}. 

According to the discussions above, one is able to generate the inphase and quadrature 
components of a bandpass channel model. Figure 2.3.3 shows the channel model is obtained 


7 see Appendix A.3 for details 



S D (f), I HM 2 In dB S D (f), | Hfc)P in dB S^f), \ HQco)| 2 in dB S D (f), | Hfc)| 2 in dB 


2.3. FACTORIZATION OF DOPPLER POWER SPECTRAL DENSITY 




Figure 2.3.1: S D (f) & H(ju)H{-ju) at (3 m = 10° 




Figure 2.3.2: Soif) & H(juj)H(—juj) at f3 m — 40° 
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by passing two white Gaussian noise process, say {wi(t)} t > o, {wQ(t)}t>o, through the rational 
transfer functions, say if/(s) and Hq(s ) which have the form of H(s). What we have, after 
filters Hi(s ) and Hq(s), are the inphase and quadrature components of the channel, which 
we denote as I(t, r) and Q(t, r). Thus, the bandpass channel model representation is 

i?e|r(t, r )e ?$ ( t,T V Wc *} = I(t, r) cos(u c t) - Q(t, r) sin(u; c t). 



Figure 2.3.3: A Bandpass Channel Model in Frequency Domain 

2.4 State Space Realizations 

In the previous section, we introduced a rational transfer function H(s ) for generating 
the inphase and quadrature components of a wireless channel, which are obtained from their 
Doppler power spectral densities. Here a stochastic state space model is used to generate 
the inphase and quadrature components. The stochastic state space model which can be 
used to realize any proper channel is given by 

x — Ax + Bu 
z = Cx + Du 
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where x G 5ft n is the state of the system, u G 5ft m is the input signal, 2 G 5ft p is the output 
signal (or say the observation of the system), A G 5ft n x 3?" is the state matrix, B G 5ft n x 5ft m 
is the input matrix, C G 5ft p x 5ft n is the output matrix and D G 5ft p x 5ft m is the direct feed 
through matrix. Here $ft denotes the set of real numbers. The block diagram is shown in 
Figure 2.4.4. 



Figure 2.4.4: A Stochastic State Space Model 

Although the state space realization for an arbitrary rational transfer function H(s ) is 
not unique, two common ways of the state space realization of a rational transfer function 
are introduced in this thesis, for generating the inphase and quadrature components of a 
wireless channel. They are so called the controllable and observable method, respectively. 
For the rational transfer function 

k 

H< - S) = s 2 + 2 Cw„s + < 

a second-order stochastic differential equation is related to the state space realization. This 
differential equation is given by 

x(t) + 2£cu n x(t) + = kw(t), x(0), x(0) given. 


where {w(t )} t > 0 is a Brownian motion. The state space realization of the Controllable 
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Canonical Form is 8 



0 1 


'o' 

X = 

-col -2Cu n 

x + 

k 


z = [ 1 0]x 


the Observable Canonical Form is 9 


x = 

2 == 

Consider a scattering environment, in which the line of sight (LOS) or specular component 
is present in the received signal. We assume that there are total N multipath components 
appearing in the received signal at time t, which one of them is the LOS corresponding to a 
Ricean channel, while the other N — 1 multipath components are considered to be Rayleigh 
channels. For each multipath component, there is a state space representation corresponding 
to the channel inphase and quadrature components. The state space wireless Rayleigh and 
Ricean fading channel models will be introduced in the following two subsections. 


0 

2 

—or 


k 

1 

n 

X + 

0 

0 

llx 




2.4.1 State Space Realization of Rayleigh Channel 

Here we will show how to construct the inphase and quadrature components of the i th 
path state space Rayleigh channel model. The rational transfer function, 

k{r) 


Hl{s) 


Le{I,Q}, 


S=JLJy 


s 2 + 2C(r)o; n (r)s + w 2 (r) | 
can be related to the following two second order differential equations, 

xi(t) + 2C(r)o;n(r)x/(t) + w 2 (r)x/(t) = fc(r)w/(t), x/(0), x/(0) given, 


x Q {t) + 2( (T)w„(r)x Q (t) + u >l(T)x Q (t) = k{r)w Q (t), xq(0), xq( 0) given, 
for constructing the Rayleigh channel inphase and quadrature. Here {wj(t)} t >o, {wQ(t)}t> o 
are two iid white Gaussian noise processes with density iV(0; 1). Towards introducing a state 
space representation, define the phase variables as follows, 

Xl x = X/, -Xqj = Xq, 


X h = x 7 , 




A . 

= XQ, 


x, = {Xj, X h }‘\ X Q = [X ai X Q ,]». 


8 see Appendix A.4 for details 
9 see Appendix A.4 for details 
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Then, the state space Controllable Canonical Form of Hj(s), Hq{s ) are 
Xj{t) = -4/(r)X/(t) + _Bj(r)tyj(t), X/(0) e 3ft 2 , 
X Q (t) = .4g(r)X 0 (() + B Q (r)w Q (t), X q (0) S K 2 , 

where 


(2.4.9) 


Ai(t) = ^q(t) = 



1 

-2C(r)a; n (T) ’ 


5/(r) - 5q(t) 


0 

/c(t) 


Let /(t) = xi(t), Q(t ) = XQ(t), then 


/(<) = [i o]X/(t), g(t) = [i o]x 0 (t). 


Define 


AiXp) 0 

B — 

B/.(n) o 

X- — 

x u (t) 

. A 

on. — 

w h {t) 

0 AqXti) 

J — 

o Bq,(t i) 

5 XV l 

X Ql (t). 

, W x 

w Qi {t) 


Then, the i th path Rayleigh channel inphase and quadrature representation of the Control¬ 
lable Canonical Form is 

Xi = AiXi + Bm Xi{ 0) 6 3ft 4 . (2.4.10) 

Thus, 


Ii(t) = [1 0 0 0 }Xi(t), Qi(t) = [ 0 0 1 OpC(t). 


2.4.2 State Space Realization of Ricean Channel 


Let {I R (t)} t > o, {Q R (t)}t> o be the inphase and quadrature components of a Ricean chan¬ 
nel. Then, the Ricean channel can be generated by assuming E[I R (t)] = r 0 cos(u 0 t + do), 
E[Q R (t)] — ro sm(<jj 0 t + 6o) (see [3]). Thus, the Ricean channel of the Controllable Canonical 
Form can be represented as follows, 


Xp(t) = Xi(t) + r 0 cos (u 0 t + do) 
Xq(t) = X Q (t ) + r 0 sin(o;ot + do) 


(2.4.11) 

(2.4.12) 


where the parameters (r 0 , w 0! 0 O } correspond to the specular component or LOS. For the j th 
multipath component, define 
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A ij(o) 0 

B — 

'B,,(t,) 0 

0 A Qj ( Tj ) 

7 

0 B Q,( t s). 


Then, the Ricean channel inphase and quadrature representation of the Controllable Canon¬ 
ical Form is 10 


Xf = AjX« + f lo 


+ BjWj, Xf( 0) € 5ft 


(2.4.13) 


where 


Thus, 


-r 0 u 0 sin(u 0 t + 9 0 ) ‘ 
r 0 ul cos(uj 0 t + 6*o) 
r 0 o;o cos^t + 9 0 ) 

r 0 u>n sin(aj 0 f + 6 0 ) . 


I*(t) = {1 0 0 0]Xf(t), Qf{t) = [0 0 1 0]Xf(t). 


Suppose the received signal over [0 T\ involves N multipath components of which one is 
the LOS. We assume that the first path is the LOS which corresponds to the Ricean channel. 
Thus, the state space representation of the output {y(t)} t >o corresponding to (2.1.2) is 

y = Y / G i X i + Dv , (2.4.14) 

t=l 

where G* = s(t — Ti)Ci(t), 

Ci(t) = [ C0 s(u; c f) 0 — sin(a ) c t) oj, D = [ cos(u c t) — sin(a; c i) j and v = 


vi (t) 


2.5 State Space Channel Model Simulations 

In this section, several simulations of Rayleigh and Ricean state space channels will be 
performed using the simulink package from Matlab. 


2.5.1 Rayleigh Channel Simulations 

- Here, a fast fading channel is simulated using model (2.4.10) as shown in the block 
diagram of Figure 2.5.5 11 . There are two possible methods. The first one is based on 
10 see Appendix A.6.1 for details 

11 Here only one of the multipath components is simulated, we choose the i th path state space Rayleigh 
fading channel model (2.4.10). It is convenient to remove the subscript i in this subsection. 
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specifying the channel parameters { u n , £, k} in prior, the second one is based on extracting 
the channel parameters from measurement data. 

Method 1: Given the values (3 m , v, f c and Eo, one can determine the state space model. 
For example, when jd m = 10°, v = 120 km/h, f c - - 900 MHz and Eq = 2, we obtain, from 
(2.3.5) to (2.3.7), w n — 637.26 rad/s, C = 0.169 s _1 , k — 2.3 x 10 4 . Figure 2.5.6 shows 
the inphase {7(t)}t>o, quadrature {Q(t )} t > 0 , phase {$(f) = tan -1 (in radian) and 

envelope (r(t) = \Jl(t) 2 + Q(t ) 2 } t >o (in dB). 

A typical Rayleigh fading envelope at 900 MHz and 120 km/h is reported in [12], 
page 173. In order to reach the same signal level, we adjust the channel model parame¬ 
ters {u> n , (, k} 12 for (3 m = 10°, v = 120 km/h, f c = 900 MHz and Eo = 2, which correspond 
to w n — 637.26 rad/s, £ = 0.169 s -1 , k = 6.11 x 10 3 (see Figure 2.5.7). Our further 
investigations, in this thesis, are based on this modified channel model. 

Method 2: Given the steady state mean of a channel envelope, namely, lim^oo E[rdB{t)\ = 
a dB, which is extracted from measurement data, one can determine the parameters {C, u> n , k} 
of the channel model as follows. Clearly, since {I(t)}t> o and {Q(i)}t>o are hd, then 

£[r 2 (t)] = E[I 2 (t) + Q 2 (t)] = E{l\t)} + E[Q\t)} = 2 E{P(t)) = 2 £[Q 2 (i)]. 

The solution of the state space model (2.4.9) is 

X L {t) = e ALt X L {0) + f e A ^-^B L w L {r)dT. 

Jo 

The mean vector is 


£[*i(t)] = e^'E^O)]. 

The covariance matrix is 

S L (t) = E[X L (t)X‘i(t)} 

= e A '- , E{X L (<S)X'[(0)]e A '^ + S‘,e A ' 1 ‘- r >B,,B , [e A ""- T) d7 l 

which satisfies 

± L {t) = A L E L {t) + E L (t)A» + B L Bi r , S L (0) - E[X L ( 0)Xf(0)], 

where 


(2.5.15) 


(2.5.16) 


12 Multiply k by 
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Rayleigh Fading Channel 



Figure 2.5.5: A Bandpass Rayleigh Fading Channel Model 






Figure 2.5.6: The Rayleigh Fading Channel; f c — 910 MHz, v = 120 km/h 
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.2 I-1-I-L-!-1 

0 50 100 150 200 250 

[a] inphase. Time [ms] 



[b] quadrature. Time [ms] 



Figure 2.5.7: The Modified Rayleigh Fading Channel; f c = 900 MHz, v = 120 km/h 


SiW = 


' E [/*(*)] 

Emm 

A 

7n(f) 

712 (f) 

E\mm) 

E[P(t )] . 


. 712 (f) 

722 (f) 


L£{I,Q}. 


The solution {E(f)} t > 0 can be implemented recursively using (2.5.16), or obtained explicitly 
through (2.5.15). 

Since a Rayleigh R.V. r(i) is generated by adding two iid zero mean Gaussian R.V.’s, 
then E[r(t)\ = a{t ). Further, if lim^^ E[r(t)] — (e.g., variance converges), then 

from the data we are given lim^oo E[r(t)} = In our case, the inphase and quadrature 

/(f), Q(t) are generated through the state space model (2.4.9). Therefore, we should find 
{k,u n ,C} so that the measured value lim^oo E[r(t)} — a = a which corresponds to 

yi Jmm- 

For a Rayleigh channel, we set E[X L (0)] = 0, e.g., X L is a zero mean Gaussian vector. 
As E[I 2 (t)} = 7 n (t) = E[Q 2 (t)}, then {fc,ai n ,(} should be chosen so that the steady state 
limt_ > oo E[I 2 (t)} = limt^oo E[Q 2 (t)} converges to a dB. Towards presenting the recipe of 
determining {k,w n , (}, we shall present the explicit solution which requires the following 
preliminary calculations. For Al corresponding to the Controllable Canonical From, we 
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have 


e ALt = ^{(sI-Al)- 1 } 

ai(t)e~^ Wnt + £*2(t)e~ < *‘ Jnt ol$ (t)e~^ nt 

Q!4 (t)e~^ Wnt ai(t)e~^ Wnt — a 2 (t)e _ ^ nt 


where s = ju , 

Qi(t) = cos (cj n Vl - C 2 t), 


a 2 (*) = sin ( Wn ^ /1 ~ C^), 

« 3 (t) = ^ 0_ c2 sin (w n Vl - C^), a 4 (t) = - sin (oj ny /1 - ( 2 a) • 


Then 


7n (*) = ( 


722 ( 0 ) + 7ii(0)£ai n 


k 2 


+ 


u, n JT= C 5 4 w 3 v ^!) Sin ( 2lJ "l/ r ~?‘)^ 2C “”‘ 

7,i(0)Cu.; + 2712(0)^+ 7 2 2 j 0 ) sin 2 

W«( 1 - C 2 ) v v ' 


A : 2 


+ 7n(0) cos 2 ( W „f^t)e- 2 W - ————ye- 2 ^" 4 


712 (A) = ( 


+ cos ( 2 ^»v / i-C*) e ' 2< ""‘ + 

722(0) - 7ii(0)cj 2 n . 


k 2 


it 2 


) sin (2„„^n^)e- 2 <“-* + 4^ (1 _ 


2n>„ v , W I 

_ 7ll(0VnC + 7l2(0)C 2 ^n + 7l2(0)^n + 722 (0)C ^2 ^ ^/l - £2^ e ~^nt 
^Tl(l _ C 2 ) 

+ 712(0) cos 2 (w„\/l - < 2 i)e“ 2< ""‘ - *2^ COS (w„/f - C 2 i)e" 2c "-‘, 


722(0 = - ^(0K±| 2 (0)C ) sin 

+ ? 1 i 0 >l + - 2 ^ 2 2 ) 0 ) + 722 . (°) - C ? . sin 2 
+ 722(0) cos 2 ( w „\/r75 t )e- 2 <-‘ + 

+ _^!_cos (2. ,/T37J t )e- 2 ^‘ + (1 - 2C* + C 2 + (2C S - 

v/T^C 1 ""“V 1 i*J e + 4Cu>„(l-C 2 ) 


As t -> 00, we can compute the steady state covariant matrix which is 




_ 


r k2 

0 i 

£ l = lim E L (i) = 

t—>00 

7n 

_7l2 

712 

722 

- 

4C w n 

. 0 

(i-2C 4 +C 2 +(2C 3 -C)V /l -C 2 )fc 2 

4Cu> n (l-C 2 ) J 



2.5. STATE SPACE CHANNEL MODEL SIMULATIONS 


27 


Since (see [19]) 


r dB = 201og 10 r = clnr 2 ; 
the mean value of the clB-record is 

roc 

E[r dB } = / clnr 2 f(r)dr, 

Jo 

where c— 10/In 10, f(r) is pdf. Then, for a Rayleigh distributed envelope, we have 

E[r dB ) = 101og 10 (2cr 2 ) - 2.51. (2.5.17) 

In our case, a = E[r dB \■ That is, 

in 0 + 2-51 , 2 

2 10 io k 

a =^- = 7ll = ic^- 

Thus, the parameters {k,w n ,(} satisfy the following equation 

k = oj n \j2C > u) n \^ 9i± ^ ■ (2.5.18) 

Figure 2.5.8 [a] corresponds to the mean of the channel envelope a = —10 dB, which implies 
k = 4.016 x 10 3 while C = 0.1691 s _1 , u n = 644.345 rad/s\ Figure 2.5.8 [b] shows the mean of 
the channel envelope when a = 0 dB, which gives k = 1.27 x 10 4 for selecting £ = 0.1691 s -1 , 
u) n = 644.345 rad/s; Figure 2.5.8 [c] displays the mean of the channel envelope at a = 10 dB , 
which means k — 4.016 x 10 4 for choosing ( = 0.1691 s _1 , u; n = 644.345 rad/s ; Figure 2.5.8 
[d] shows the mean of the channel envelope at a = 0 dB, which gives k — 5.34 x 10 3 while 
( = 0.1 a" 1 , u> n — 250 rad/s. 

2.5.2 Ricean Channel Simulations 

This is very similar to the Rayleigh case, because the Rayleigh case is related to Ricean 
by (2.4.11) and (2.4.12). A flat fast fading channel is simulated using the model (2.4.13) 
as shown in the block diagram of Figure 2.5.9 13 . There are two ways of Ricean channel 
simulation. Both of them are based on a prior knowledge of the LOS parameters {r 0 , u/ 0 ,6b). 

Method 1. Quite like the Method 1 of Rayleigh. Given LOS values {r 0 , ui 0 , $ 0 } and (3 m , v, 
f c , Eq, one can determine the Ricean state space channel model. For example, when j3 m = 


13 Here again only one multipath component is simulated. We select the j th path (j = 1) Ricean channel 
model (2.4.13). It is convenient to remove the subscript j in this subsection. 
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Figure 2.5.8: The Envelopes of Rayleigh Fading Channel; f c = 910 MHz, v — 120 km/h. 



Remark: A_l = A_Q, B_l = B_Q, C_l = C_Q 


Figure 2.5.9: A Bandpass Ricean Fading Channel Model 
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10°, v = 120 km/h, f c — 910 MHz and E 0 = 2, we obtain w n = 644.34 rad/s, ( = 0.169 s 
k — 6.21 x 10 3 . Figure 2.5.10 shows the inphase {/ J? (t)}t>o ) quadrature }«>o, phase 

{$(£) = tan" 1 (in radian) and envelope (r(t) = y / I R (t) 2 + Q R (t) 2 } t > 0 (in dB) 

for ro = 1, loq = 6°, do = 30°. 

Method 2. The steady state mean of the square envelope of a Ricean channel is 

£[r 2 ] = 2a(0) + r 2 , 

which we introduced in the chapter 1. Given the steady state mean of a Rayleigh channel 
envelope in dB, say 


lim E[r dB (t )] = a, 

t—>oo 

which is extracted from measurement data, then by using (2.5.17), we obtain 

E[r 2 ] = 10° + io +^o- 

We already know that given the steady state mean of a Rayleigh channel envelope, one can 
determine the state space Rayleigh channel model parameters {C, w n , k} from (2.5.18). From 
the above equation, one can determine the steady state mean of the square envelope of the 
Ricean channel by given the LOS parameters {r 0 , Uq, do}- For example, while r 0 = 1, Uq = 6°, 
d 0 — 30°, Figure 2.5.11 [a] corresponds to the mean of the square envelope is E[r 2 ] = 1.42 dB 
for a = — 10 dB, which implies k = 4.016 x 10 3 while C — 0.1691 s" 1 , u n = 644.345 rad/s ; 
Figure 2.5.11 [b] shows the mean of the square envelope is E[r 2 ] = 8.88 dB when a = 0 dB, 
which gives k = 1.27 x 10 4 for selecting ( = 0.1691 s" 1 , ui n — 644.345 rad/s] Figure 2.5.12 
[a] displays the mean of the square envelope is E[r 2 ] = 25.49 dB at a = 10 dB, which means 
k — 4.016 x 10 4 for choosing £ — 0.1691 s" 1 , u n = 644.345 rad/s; Figure 2.5.12 [b] shows the 
mean of the square envelope is E[r 2 ] — 8.88 dB when a = 0 dB, which gives k = 5.34 x 10 3 
while C = 0.1 s" 1 , u) n = 250 rad/s. 
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Figure 2.5.10: The Ricean Fading Channel; f c = 910 MHz, v = 120 km/h 
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Figure 2.5.11: Part 1: Ricean Channel Square Envelopes; f c = 910 MHz, v = 120 km/h 




[b] a =0 dB. Time [ms] 


Figure 2.5.12: Part 2: Ricean Channel Square Envelopes; f c — 910 MHz, v = 120 km/h 
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Chapter 3 

Channel Model Estimation 


These results of chapter 2 are first employed to construct wireless state space Rayleigh 
and Ricean fading channel models. The objective of this chapter is to employ least-square 
estimation techniques to estimate the channel models from noisy received signals. The focus 
on is based on Kalman filtering theory, to extract the channel inphase, quadrature and square 
envelope estimations. The following items will be discussed 

1. Continuous time channel state estimation. 

2. Discrete time channel state estimation. 

3. Computation of the conditional Rayleigh envelope and phase distributions. 

All investigations will be done by assuming that the parameters {A, B, C, D} of the state 
space dynamic system (2.4.8) are known. If {A, B, C, D} are unknown, one can also employ 
the state space estimations to estimate the parameters {A, B,C, D} as in [7] and [9]. This 
technique is so called systems identification. 

3.1 States Estimation 

A general state space dynamic system and measurement models are given by the following 
equations 

X = AX + Bw 
y — CX + v 

Here A'(O) ~ iV(/x, E), w ~ N( 0, Q ), v ~ N( 0, R ), and {u;(f)} t >o, {v{t)}t>o are independent 
white Gaussian processes, which are uncorrelated with -X’(O). The algorithm corresponding 
to the Kalman filter is given by the following equations. 
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Initialization: 


P(0) - Var(X( 0)) = E(0), X(0) = X(0) = Mean(X( 0)) = /q 

Error Covariance (Riccati Equation): 

P(t) = AP(£) + P(t)A tr + BQB tr - P(t)C , ‘ r P“ 1 CP(0; 


Filter Gain: 


K(t ) = P(i)C* r P -1 ; 


Estimation (Filter): 


X = AX + K(t)(y-CX). 


3.1.1 Continuous Time 


First we apply the Kalman filter algorithm for a flat fading channel and then we introduce 
the multipath Kalman filter algorithm. The Kalman filter corresponding to (2.4.10), (2.4.13) 
and (2.4.14) is given below. 

Rayleigh Channel Kalman Filter. On the i th path (i ^ 1), the state space dynamic 
system and measurement models are given by (2.4.10), (2.4.14), respectively. For simplifying 
the notations, we represent the i th path state space Rayleigh channel model as follows, 


where 


X = AX + Bw 
y = GXP Dv 


(3.1.1) 


X=X„ A = A,, B=B h w=m, G = Gi, 


and Dv is a white Gaussian noise added on the i th path. Since v has a coefficient D, we should 
substitute P” 1 by DR~ l D tr in the Kalman filter equations above 1 . Thus, the least-square 
estimates of {I(t)}t> o, {Q(£)}t>o and {r 2 (£)}t>o are given by the following equations. 
Initialization: 


P( 0) = Var(X(0)), X(0) = X(0) = Mean(X{ 0)); 


1 see Appendix A.5.1 for details 
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Error Covariance (Riccati Equation): 

P(t) = AP{t) + P{t)A tr + BQB tr - P(t)G tr (t)(D(t)i? J D tr (t))- 1 G(t)P(t); (3.1.2) 
Filter Gain: 


K(t) = P{t)G tr (t)(D{t)RD tr (t))-\ 


where 


Estimate (filter): 


R± 


<T 2 0 

A 

1 

o' 

0 a 2 

and Q — 

0 

1 


X = AX+ K(t){y(t)-G(t)X ); 
Inphase and Quadrature Estimates: 


(3.1.3) 


m = e 


Q(t) = E 


I{t)\{y(s)] 0 < s <t} 


Q(t)\{y(s)-, o <s<t} 


= 10 0 o }X(t), 


= [001 0]X(t); 


Square-Envelope Estimate: 


r 2 (t) = P(t ) + Q 2 (t ) + e]{t) + e 2 Q (t), 

where e]{t) = £[(/(£) - J(£)) 2 ], e^(t) = i?[(Q(t) - Q(t)) 2 ] are the mean-square errors, which 
are obtained from (3.1.2), that is, the elements Pn{t), ^ 33 (t) of P(t), respectively. 

Ricean Channel Kalman Filter. On the j th path (j = 1), the state space dynamic 
system and measurement models are given by (2.4.13), (2.4.14), respectively. Similar to the 
Rayleigh, the j th path state space Ricean channel model is 


X — AX + fi os + Bw 
y = GXP Dv 


(3.1.4) 


where 


X±X; 


R 


A A 7—) A j-j . A . .—y 

A = Aj, B = Bj, w = Wj, G = Gj, 


and Dv is a white Gaussian noise added on the j th path. The least-square estimates of 
{I R (t)} t > 0 , {Q R {t)}t> 0 and {r^(t)} t >o are given by the following equations. 
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Initialization: 


P{ 0) = Var(X(0)), X(0) = X(0) = Mean{X{ 0)); 

Error Covariance (Riccati Equation): 

P(t) = AP(t) + P(t)A tr + BQB tr - P(t)G ir (t)(L>(t)i? J D tr (t))- 1 G(t)P(t); (3.1.5) 


Filter Gain: 


K(t) = P(t)G tr (t){D{t)RD tr (t))-\ 


where 


R = \“ v 


al 0 


0 ai 


and Q = 


1 0 
0 1 


Estimate (filter) 2 : 

t = AX + f los + K(t)(y{t) - G(t)X); 
Inphase and Quadrature Estimates: 


(3.1.6) 


I R {t) = E 


Q R {t) = E 
Square-Envelope Estimate: 


I R {t)\{y{s); 0 <s<t} 


Q (t)\{y(s)-, 0 <s <t} 


1 0 0 0 }X(t), 


= 001 0 ]X{t)] 


r 2 R {t) = + Q R2 {t) + Cjji(t) + e 2 Q R (t), 

where e] R {t ) = E[(I R (t)-I R (t)) 2 ], e 2 QR {t) = E[(Q R (t)-Q R (t)) 2 } are the mean-square errors, 
which are obtained from (3.1.5), that is, the elements Pn{t), ^(i) of P{t), respectively. 

Multipath Channel Kalman Filter. Consider a multipath channel having N multi- 
path components, the first of which is a Ricean channel, while the remaining N — 1 com¬ 
ponents are Rayleigh channels. This multipath channel model can be represented by using 
(2.4.10), (2.4.13) and (2.4.14), as follows. 


2 see Appendix A.6.2 for details 


X = 
y = 


AX + Fi os A Bw 
GX + Dv 


(3.1.7) 
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where X(0) G K 4iV , 

' A x ■ ■ ■ 0 1 \ B l •■• O' 

a= : ; , g = [g x 

0 ■ • • An 0 ■ • • Bn 

and 

Fi0S = L// 0 r s 0 0^0 o' 0 0^0 oJ ' 

Thus, one can obtain the least-square estimates of the inphases, quadratures and square 
envelopes by the following equations. 

Initialization: 

P{ 0) = Var(X( 0)), X(0) = X(0) = Mean{X{ 0)); 

Error Covariance (Riccati Equation): 

P(t) = AP(t) + P(t)A tT + BQB tr - P(t)G tr (DRD tr )~ l GP(t)] (3.1.8) 

Filter Gain: 

K(t ) = P(t)G tr {DRD tr )~ l , 

where 

'<%■■■ oi r i • • • o■ 

R = : : , Q= : , R, Q G K 2iV x U 2N - 

0 • • • a 2 v \ [ 0 • • • 1 

Estimate (filter): 

t = AX + Fi os + K(t)(y(t) - GX), (3.1.9) 

Inphase and Quadrature Estimates: 

%{t) = 

i 

= 0 0 0 0 ••• 'l o" 0 o' •••0 0 0 0 ] x ( t ), 

Q,(t) = E Qi(i)l{»(«)}i=o 

i 

= 0 0 0 0 ••• 0 o" 1 o' ••• 0 0 0 0 ] x ( t ). 
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Here I\(t) denotes the estimate of I R (t), Qi(t) denotes the estimate of Q R (t). 
Square-Envelope Estimate: 

rKt) = ^(t) + QKt ) + 4(*) + e 2 Qi (t), 

where e%(t) = E[(Ii(t) — Ii(t)) 2 ], e^f) = E[(Qi(t) — Qi(t)) 2 ] are the mean-square errors, 
which are obtained from (3.1.8), that is, elements of PAi-sAi-zit), PAi-iM-i(t) of P(t), 1 < 
i < N, respectively. 


3.1.2 Discrete Time 


Discrete Time Multipath Channel Kalman Filter. Although estimations can be 
done by using continuous time model, we prefer the discrete time model in this thesis. In 
order to use the discrete time algorithm of the Kalman filter, we need to do the discretization 
for the continuous time multipath state space channel model (3.1.7). Mathematically, we can 
not prove the existence of a white noise, but we can represent a white noise in terms of the 
velocity of a Brownian motion, that is, v(t) — dw(t)/dt = w(t). Along this way, we can obtain 
a very good result after the discretization. Sampling (3.1.7) on a time interval [a b] by using 
the Ito process method with respect to an equal time subinterval [t n t n+ 1 ], Vf n £ [o b ], 
then we obtain the discrete time multipath state space channel model as follows 3 

X n+ i = rA n + T + AAw n 
y n — LlX n + ^Lv n 

where 


(3.1.10) 


T = 5A + I, T = SF los , A = B, 


n = 


S(t n - T^C^tn) 


0 


0 


s(t n - r N )C N (t n ) 


vi(t n ) 
VQ^tn) 


^ = [cos(w c f n ) — sin(o; c t n ) ], v n 
Here 5 = t n+ \ — t n is called the step size, I G 3ft 4iV is the unit matrix, 


Ci(t n ) = [cos(u c t n ) 0 -sin (uj c t n ) 0] 


3 see Appendix A.5.2 for details 
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and A w n = w(t n+ 1 ) — w(t n ) E 9? 2JV are iid white Gaussian noises, which each element 
of Aw n ~ N(0;5). Thus, the algorithm of the discrete time Klaman filter produces the 
least-square estimations of the inphases, quadratures and square envelopes are given by the 
following equations. 

Initialization: 

P 0 = Var(X o), X 0 = A 0 = Mean( X 0 ); 

Error Covariance (Riccati Equation): 

p n+ 1 = rp n r 4r + AQA tr - rp n Q tr (rtp n o tr + *ira* r )- 1 fip„r tr , (3.1.11) 

where 

' <5 O' 

Q= : j E^ N ] 

° ••• 5 _ 

Filter Gain: 

K n - r P n O tr {OP n n tr + 

Estimate (filter): 

X a+ i =TX n + r + K n (y n - SlX n ); (3.1.12) 

Inphase and Quadrature Estimates: 

ii, = E liJfoUto 

i 

= [o 0 0 0 ••• \ 0^0 o' •••0 0 0 o]x n , 

Qin ~ E Qin\{ykYk -0 

i 

= 0 0 0 0 ••• 0 o"i o' •••0 0 0 o ] x n . 

Here I\ n denotes the estimate of Ifi, Qi„ denotes the estimate of Q^. 

Square-Envelope Estimate: 


rl-lf n + Ql + ej in + el, 



40 


CHAPTER 3. CHANNEL MODEL ESTIMATION 


where ej in = E[(I in - J in ) 2 ], e 2 Q in = E[(Q in - Qi n ) 2 \ are the mean-square errors, which are 
obtained from (3.1.11), that is, elements of -P( 4 i- 3 , 4 i_ 3 ) n , P{a-i,u-i) n of P n , 1 < i < N, 
respectively. 

Similarly, we can find the discrete time state space Rayleigh channel model of the i th 
path and Ricean channel model of the j th path using the Ito method. 

3.2 Conditional Rayleigh Envelope and Phase Distri¬ 
butions 


Channel envelope and phase are two very important parameters in channel estimation. 
Since the conditional densities of the envelope and phase can be used to compute any condi¬ 
tional estimation of the functions of r(t ) and 9(t), we should obtained their explicit expres¬ 
sions. We are interested in computing the conditional distributions of the channel envelope 
and phase, given the noisy data. Least-square estimation of Gaussian processes state that 
the optimal estimates of Xj{t ) 6 3ft 2 , Xq(t) e 3ft 2 given (y(s); 0 < s < t} are the conditional 
means 

% (t) = E{X,(t)\YS\, % (i)e» 2 . 

* 3 (t) = E\X Q (t)\Y’l X 5 (f)e 3 t 2 . 


where Yq = (y(s); 0 < s < t}. The conditional joint density of the inphase and quadrature 
is still Gaussian 4 , given by 

f(x h ,x Ql ,t\Y*) = dx ® Qx ^ Prob(X h (t) < x h ,X Ql (t) < x Ql |l£) 

where 

*=[*'■ l = 

_ X Q i _ 

{M(t )} t >o is the conditional covariant matrix of {x(t)} t > o, and each element of M(t) is 
related to the Riccati equation (3.1.2) as follows, 


X,,(t) 


E\mw] ' 

x Qt (t) 


E[C(()|r 0 ‘] 


M(t) = 


Pu(t) 

Pn(t) 


PUt) 
p 33 (t) ’ 


Pis(t) = P,i(t). 


4 see Appendix A.5.3 for details 
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3.2.1 Conditional Envelope Distribution 


Let Z — X] x + Xq x ; then we have, by definition, 

Prob(Z(t) < z\Y*) = Prob(Xl(t) + X 2 Ql (t) < z\Y$) 

= L , f{xi 1 ,x Ql ,t\Yo)dx h dx Ql 


■L 


x^+x^Kz 


x\+x 2 Ql <z 2irJ\M(t 




Define the inverse matrix of M(t) 


M~\t) = 


p33(t) P 13 CO 


Pn(t)PM - P^t) [-Pvi{t) Pn(t) 
A [ s(0 m(i) 1 


m(t) n(t) ’ 


and introduce to the spherical coordinates 


x h = y/pcosd 
x Qi = v^ sin0 


Then 5 


n 27r 1 

-/(p,M^o)<W0 


r _i 

jo o-tt./h 


(t)-m(t)Xj 1 (t)XQ 1 (t)-^-p 


2xJ M t) 


oo / _1\ko —j h 

X £ r (9(t)X h (t) + m(t)X Ql (t)) 

h,i,j,k=0 J 

x(n(t)XQ l (t) + m(t)X h (ii'j : (gU) - n(t)ym k (t) 

»^ r(Mf±i)r(!±2at±i) 

^ p( M-t+ 2 , 7 + 2 fc+ 2 ) u f J ‘ 

Thus, the conditional probability density function of the Rayleigh channel envelope is 


fR(r(t)\YS) 


0 -‘Mx Jl \t)-^X Ql \t)-mWi l (.t)XQ 1 (t)- 3 % 1 r 2 


7T \/\M(t) 


■ ( 1 p/ M-fc +1 t+frf+fc-fl N _ 


2 V 2 
< h-|-i-h2j-h2fc-|"2 > 




5 see Appendix A.5.3 for details 
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3.2.2 Conditional Phase Distribution 

Here we should first introduce an expression for the conditional density of tan - i m) 
given { y(s ); 0 < s < t}. Define 

{ x^ = p cos 9 
x Ql = psin# 

and then 9 = tan -1 Hence 6 

Prob(@ < <P\Yq) = Prob^n- 1 ^ 

x h 

= / . \ f( x h > x Qi > t\Y^)dx h dx Ql 

= /°° [* - -L== e~^ (t)) tr M~ 1 (t) (x-xJ^q (*))pdQdp 

Jo J o 2n s J\M(t)\ 

f<t> 1 
■1° 2irJ\M(t)\ 


-i&X h \t)~^X Ql \t)-mX h (t)X Qx {t) R\B) 


X 


J_ p -R 2 (B) + 


2 S\9) 2 S 2 (0) 


')dd, 


where 


S(6) = 


l9(t) g(t)-n{t) . 2 


sin 2 0 + m(i) sin 0 cos 9, 


m = 


(g(t)X h (t) + m(t)XQ 1 (t)') cos 9 + (n(t)XQ x (t) + mi^Xj^t)^ sin 9 


2 S{9) 


Thus, the conditional density of 9(t) = tan -1 given (y(s); 0 < s < £} is 
f(9(t)\Y t ) = 1 + R{0)y/ne R {e) r ^lE;\t)-^x^\t)-rn(t)xP(t)x^(t) 


6 see Appendix A.5.3 for details 



Chapter 4 

State Space Channel Estimation: A 
Case Study 


In this chapter, we present a case study of channel state estimation for flat and frequency- 
selective channels. The estimation approach is based on the Kalman filter presented in 
the previous chapter. The following items, inphase, quadrature and square envelope, are 
estimated. For the flat fading channels, we investigate both of slow and fast fading cases. 
For frequency-selective channels, only the slow fading case is investigated. 


4.1 Frequency-Nonselect ive Fading Channels 

A frequency-nonselective fading channel is also called a flat fading channel. In this 
section, one Ricean channel model given by (3.1.4) and one Rayleigh channel model given 
by (3.1.1) are performed, respectively. In order to address the performance of the estimator, 
the received signal-to-noise power is computed using the corresponding state space model. 


4.1.1 Rayleigh Channel 


Let X(0) ~ N( 0; E(0)), then, theoretically, the state space channel model (3.1.1) has the 
property of Rayleigh fading nature. 

(i) Rayleigh Channel Received Signal-to-Noise Ratio. The received SNR is 
defined by 


SNR = 


k E 

1 G(t)X(t) I 2 dt 

y e 

Jo* 1 D(t)v(t) I 2 dt 


(4.1.1) 
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where T s is the signaling interval. There are two possible ways of using the equation (4.1.1) 
for the simulation purposes. The first method is specify the desired SNR, and then determine 
how much the average power of the transmitted signal should be. The second method is to 
specify the desired SNR, the average noise power, and then determine the power of the 
receiver. Let E s be the average noiseless received signal power. Then, 


E, = Le 


T - J a ‘\G{t)X(t) | 2 it 

= Y [‘ s 2 (t)CV)i5[x(t)V’'(t)] C* r (i)df. 


Narrow Band: If {s(t )} t >o is a narrow band, and hence without loss of generality is assumed 
to be s(i) = 1, then by performing the above integration, we get 


E x = 


l r T <- 


J o 3 C(t)E[X(t)X tr {t)]C tr (t)dt 
_ 2 C 2 A - 2p, -fa- @2 2Cu n T* , Q 

4C WT. 


+ 


+ 


+ 


2 CAv / l ~~ C 2 + 2 A + 2C 2 & + ~ ( 2 p2 + (3 3 

4(,L0 n T s 


Cft - C/?3 - 2C/?s - 2p ly /T-^ 
4 u n T s 


cos {2uj n ^l - C 2 T s )e 2 ^ UJnl 


2 V T^^ + x /T^/l3 + v / r : C ? /l2 . 


4 u) n T s 

where 7 ( 0 ) = E jX(0)X tr (0)] and 

'722(0) +7n(0)Cw n k 2 


a-(- 


a>nV 1 - C 2 

Pz = 7n(0), /?4 = 


k 2 


4Cu;3(l-C 2 )’ 
Then the additive noises are given by 1 


), & = 

/% = 


sin (2u n \jl - C 2 T s )e 2CunTs , 


7n(0)C 2a; n + 27 i 2 ( 0 )C^n + 722 ( 0 ) 

^(1~C 2 ) 

a 2 a _ fc 2 

> Pe A/ >.. 3 • 


4wJ(l - C 2 ) 


4C^ 


vi ” 


E s 

m 

V Q. 

\ 

, n ( SNR) dR 

10 10 

n Q 


where {rr/(t)}t>o ~ AT(0,1), {r< 3 (t)}t>o N( 0,1) are independent white Gaussian noise 
processes, and the covariance matrix is given by 


1 see Appendix A.7 for details 
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1 0 
0 1 

Thus, the average additive noises power is 


Rv 


E s 


(SNR) dB 

10 10 


1 f ft 

En=^E / 
T a |yo 


D(t)v(t) j 2 dt 


= i ft D(t)EUt)v*(t)\D‘‘'(t)dt = 

J 3 JO 


E. 


(SNR) dB 
10 10 


(ii) Flat Slow Fading Channel Estimation. Let f c = 910 MHz , v = 60 km/h , which 
implies a maximum Doppler shift f d — 50.556 Hz, Doppler spread Bd = 2 f d = 101.112 Hz 
and coherence time T c = 9.89 ms. Select T s — 600 [is, then we have T s « T c , and therefore 
the channel is slow fading. The results of the inphase, quadrature, square envelope estimation 
are shown on Figure 4.1.1 for ( SNR) d B = —3 dB. 
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Figure 4.1.1: Flat Slow Rayleigh Fading Channel; ( SNR) d B — — 3 dB, v — 60 km/h 


• ( SNR) d B = —3 dB, v — 60 km/h, T a = 600 [is, T c = 9.89 ms 
Figure 4.1.lja] shows the inphase and its least-square estimate. Clearly, tracking is 
achieved within the first 600 [is; Figure 4.1.1[b] shows the quadrature and its estimate. 
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Obviously, tracking is also achieved within the first 600 ^s; Figure 4.1.1[c] shows the 
mean square errors of the inphase and quadrature in dB] Figure 4.1.1 [dj shows the 
square envelope r 2 (t) and its estimate r 2 (t) in dB. Definitely, tracking is achieved 
within the first 600 fis with error. 

(iii) Flat Fast Fading Channel Estimation. Let f c = 910 MHz, v > 60 km/h, then 
T c < 9.89 ms. Choose T s = 250 ms, which implies T s » T c , and thus the channel is fast 
fading. The results of the inphase, quadrature, square envelope estimations are shown as 
follows. 

• ( SNR)dB = — 3 dB, v — 60 km/h, T s = 250 ms, T c — 9.89 ms 

Figure 4.1.2(a),(b) displays the inphase, quadrature and their least-square estimate, 
respectively, while Figure 4.1.2(c) displays their mean-square errors. Figure 4.1.2(a),(b) 
show that tracking of {I(t),Q(t)} t >o is possible with mean-square error of 0.049 for 
the inphase and 0.052 for the quadrature at 4 ms. Figure 4.1.2(d) display the square- 
envelope and its estimator. Unfortunately, tacking is not possible when the channel 
experiences deep fades at -3 dB. Next we provide simulation for 10 dB. 

• ( SNR)dB = 10 dB, v = 60 km/h, T a = 250 ms, T c = 9.89 ms 

Figure 4.1.3[a],[b] shows the inphase, quadrature and their least-square estimate, re¬ 
spectively, while Figure 4.1.3[c] shows their mean-square errors. Clearly, these Figures 
illustrate the perfect tracking properties of the estimator within 3 ms, respectively, 
with a mean-square error of —39.1 dB for the inphase, and —33.3 dB for the quadra¬ 
ture. Figure 4.1.3[d] shows the square envelope and its estimate. It is obvious that 
tracking is possible at 10 dB even for fast fading channels. 

• ( SNR)dB — — 3 dB, v — 100 km/h, T s — 250 ms, T c = 5.9 ms 

Figure 4.1.4(a),(b) show the inphase, quadrature and their least-square estimates. 
Clearly, tracking is achieved but not well because of the low received signal to noise ra¬ 
tio and the high velocity; Figure 4.1.4(c) shows the mean-square errors; Figure 4.1.4(d) 
shows the square envelope r 2 (t) and its estimate r 2 (t). The estimator is not able to 
... track the square envelope because of the low received signal-to-noise ratio and the large 
velocity. 

• ( SNR)dB — 10 dB, v — 100 km/h, T a — 250 ms, T c = 5.9 ms 

Figure 4.1.5(a),(b) show the inphase, quadrature and their least-square estimates. 
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Clearly, very good tracking is achieved, because of the high received signal-to-noise 
ratio; Figure 4.1.5(c) shows the mean square errors; Figure 4.1.5(d) shows the square 
envelope r 2 (t ) and its estimate r 2 (t). The estimator can track the fading even at points 
which of deep fading because of the high received signal-to-noise ratio. 


However, better tracking of the deep fade is expected if one employs ( SNR)dB = 15 — 20 dB. 

4.1.2 Ricean Channel 


If X(0) ~ jV(0;£(0)), then, theoretically, the state space channel model (3.1.4) has the 
property of Ricean fading nature. 

(i) Ricean Channel Received Signal-to-Noise Ratio. The received SNR for a 
Ricean channel can be obtained from that of a Rayleigh channel by using (2.4.11), (2.4.12). 
Therefore, the received SNR is defined by 


SNR r = 


y s e 

Jc?‘\G(t)X(t) + G(t)h.(t) | 2 dt 

y.e 

Jo* 1 D(t)v(t) | 2 dt 



where 


h s (t) 


r 0 cos (u 0 t + 6 0 ) 
0 


r 0 sm(u> 0 t + 9 0 ) 
0 


Similar to the Rayleigh channel, the average noiseless received signal power is 

I G(t)X(t) + G(t)K(t) P dt 

= ^ jT' S 2 (t)C(t)B[x(t)X ,r (t)]c , '-(t) [ it + i £ s 2 (t)C(t)E[h,(t)h‘;(t)]C ,r (t)dt 
= parti (t) + part 2 ( t ), 


where 

part\{t) 


l r T * 
T s 


s 2 (t)C(t)E[X(t)X tr (t)]C tr (t)dt = ^ s 2 (t)-/ u (t)dt = E s , 

1 r T s r , s 2 (t)r% s 2 (t)r„ sin (2(iv c + coo)T s + 

P ar h(t) = -/ o 5 (t)C(t)F;[fi s (t)firw]^(t)dt = -Y- +-4K+0*)r;- 


2 
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Thus, the additive noises are given by 


Vi 


E? 

ni 

V Q 

\ 


n Q 


the covariance matrix is given by 


Ef 


jdK __ 

v ~ JSNR R ) dn 


W 


10 


0 

1 


and the average additive noise power is given by 


E *=k E f 


D(t)v(t ) | 2 dt 


i £ D(t)E[v(t)v t '(t)]D tr (t)dt = 


E? 


( SNR R ) dR 

10 10 


(ii) Flat Slow Fading Channel Estimation. Let / c = 910 MHz, v = 60 km/h, which 
implies f d = 50.556 Hz, B d — 2 f d = 101.112 Hz and T c = 9.89 ms. Select T s — 600 /xs, 
then we have T s «T C , and therefore the channel is slow fading. The results of the inphase, 
quadrature, square envelope estimation are shown on Figure 4.1.6 for ro = 1, u>o = 5°, 
0o = 30° and ( SNR) dB = -3 dB. 
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Figure 4.1.6: Flat Slow Ricean Fading Channel; ( SNR) d B = —3 dB, v = 60 km/h 
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• ( SNR)dB = —3 dB, v = 60 km/h, T s — 600 /is, T c = 9.89 ms, r 0 = 1 

Figure 4.1.6[a] shows the inphase and its least-square estimate. Clearly, tracking is 
achieved within the first 600 /is; Figure 4.1.6[b] shows the quadrature and its estimate. 
Obviously, tracking is also achieved within the first 600 fis; Figure 4.1.6[c] shows the 
mean-square errors of the inphase and quadrature in dB ; Figure 4.1.6[d] shows the 
square envelope r 2 (t) and its estimate r 2 (t) in dB. Definitely, tracking is achieved 
within the first 600 fis. 

(iii) Flat Fast Fading Channel Estimation. Let f c = 910 MHz, v = 60 km/h, then 
T c = 9.89 ms. Choose T s = 250 ms, which implies T s » T c , and thus the channel is fast 
fading. The results of the inphase, quadrature, square envelope estimation are shown on 
Figure 4.1.7 for r 0 = 1, u> 0 = 5°, 6 0 = 30° and ( SNR)dB — 10 dB. 

• ( SNR)dB = 10 dB, v = 60 km/h , T s — 250 ms, T c = 9.89 ms, ro = 1 

Figure 4.1.7(a), (b) shows the inphase, quadrature and their least-square estimate, 
respectively. Clearly, trackings are achieved at ( SNR)<ib = 10 dB. While Figure 
4.1.7(c) shows their mean square errors. Figure 4.1.7(d) shows the square envelope 
and its estimate. It is obvious that tracking is possible at ( SNR)d B = 10 dB even for 
fast Ricean fading channels. 

4.2 Frequency-Selective Fading Channels 

If we transmit a narrow pulse, thus a wideband signal, through a multipath fading chan¬ 
nel, then the channel is frequency-selective in this case. The received signal model is given 
by (2.1.2), which consists of N multipath state space fading channel models given by (3.1.7). 
In this section, we employ the Kalman filter algorithm to produce the estimates of the in- 
phase, quadrature and square envelope for each multipath component. First, we show how 
to determine the multipath components N, then we simulate the frequency-selective slow 
fading channel. 

(i) Determination of the Multipath Components N. The multipath components 
N can be obtained by the formula N = [T m B s ] + 1 with the deterministic time delay {r* = 
iTs}iLi 2 - The typical values of multipath delay spread T m are listed in the table of Typical 


2 the formula can be found in [17], page 797 
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Measured Values of RMS Delay Spread, which is reported in [ 12 ], page 162. Based on that 
table, we are able to compute the multipath components N. We choose T m — 2 ps which 
corresponds to the New York city environment, carrier frequency is 910 MHz. Let the 
bandwidth of the transmitted signal be B a = 2.5 MHz, then the symbol duration (signaling 
interval) T s = 1 /B s — 0.4 ps. Thus, we get the multipath components N — 6 , the time 
delay {r; = *T s }f =1 . 

(ii) Received Signal-to-Noise Ratio. The received SNR is defined by 


SNR A 

±E[g- | D(t)v(t) I 2 dt] ' 

By simple computing, we obtain the average noiseless power of the received signal 3 


E a 


-E 


[ m I E I 2 dt] = [ m 5Z s2 (* ~ Ti )(7n + 4i)dt, 

Jo i—i T m Jo i=1 


where 7 ^ = parti(t) is the steady state variance of the inphase (quadrature) of the i th 
multipath component, d\ x — part 2 (t) (see [13]), d\ x = 0 for i ^ 1 and 


s{i - Ti) = 


const ^ <t <Ti + T s 
0 otherwise 


and the average additive noise power 


E v — — 

Trr 


E 


/ Tm | D(t)v(t) | 2 dt] = S f" D(t)E[v(t)v ,r (t)]D«{t)dt = al 
Jo 1 yn Jo 


(iii) Frequency-Selective Slow Fading Channel Estimation. The estimations are 
done by using Simulink blocks, which is referred to (3.1.7), as shown in Figure 4.2.8. This 
multipath fading channel model consists of one Ricean and five Rayleigh’s. 

Next, we are going to simulate the estimators of the inphase, quadrature and square 
envelope for each path. Each path at the receiver corresponds to the same transmitted 
signal with a certain transport delay (T.D.). Let f c = 910 MHz, v = 60 km/h, then 
T c = 9.89 ms. Thus, T s < T m < T c . The multipath fading channel shows the property 
of slow fading nature. Suppose the transmitted signal is a pulse with duty cycle 25, which 
the pulse width is 0.1 ps. Then the received signal has total 6 pulses separated by 1 /B s 
during T m . The results of the simulations are shown from Figure 4.2.9 to Figure 4.2.17 for 
(SNR)d B - 20 dB. 


3 see Appendix A.7.5 for details 
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Multipath Rayleigh / Ricean Fading Channel 



Figure 4.2.8: Multipath Fading Channel Estimator 


Figure 4.2.9[a] shows the transmitted signal; Figure 4.2.9[b] displays the absolute value 
of the received noiseless signal \y{t)\ = | ]T-Li (li(t) cos(oj c t) — Qi(t)sm(u c t)^s(t — T;)|, the 
multipath delay profile, which illustrates the amplitude of the received signal of each path 
goes down because of the longer distance and larger time delay. Path #1 is a Ricean 
channel, the steady state mean of the channel square envelop is 8.88 dB which corresponds 
to the channel model parameters k = 2.197 x 10 3 , w n = 322.1724 rad/s, £ = 0.1691 s -1 
and the LOS component parameters r 0 = 1, u) 0 = 31.4 Hz, 6 0 = 30°. Path #2 is a 
Rayleigh channel, the steady state mean of the channel envelop is 10 dB which corresponds 
to the channel model parameters k — 6.748 x 10 3 , w n = 322.1724 rad/s, £ = 0.1691 s -1 . 
The steady state means of the channel envelops from path #3 to path #6 are 15 dB, 
20 dB, 25 dB and 30 dB, which correspond to the channel models parameters ({k,u> n , £}) 
{1.2xl0 4 ,322.1724,0.1691}, {2.1339x 10 4 ,322.1724,0.1691}, (3.7947x 10 4 ,322.1724,0.1691} 
and {6.748 x 10 4 ,322.1724,0.1691}; Figure 4.2.9[c] shows the absolute value of the received 
noisy signal; Figure 4.2.9[d] displays \y('t)\ and its estimate \y(t)\ = | Z!i=i (Ji(t) cos(u c t) — 
Qi{t ) sin(w c i))s(t -r f )|. 

Figure 4.2.10 shows the received signal of path #1, path #2, path #3 and their estimates. 
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Figure 4.2.9: [a]: the transmitted signal; [b]: the multipath delay profile | y(t) |; [c]: the 
received signal with noise; [d]: | y(t) | and its estimate 

Figure 4.2.11 shows the received signal of path #4, path #5, path #6 and their estimates. 
Figure 4.2.10 and Figure 4.2.11 show that tracking is not achieved at the beginning of each 
path, but it is achieved as the time elapses. 

Figure 4.2.12[a] shows the inphase and its least-square estimate. The initial value of the 
inphase is r 0 cos(6 0 ). Obviously, tracking is achieved from the time 0.4 ys to 0.5 ys, because 
there is a transmitted signal on [0.4 0.5]; Figure 4.2.12[b] displays the quadrature and its 
least-square estimate. The initial value of the quadrature is rosin(0 o ). Clearly, tracking is 
achieved between the time 0.4 ys and 0.5 ys; Figure 4.2.12[c] shows the mean square errors 
of the inphase and quadrature; Figure 4.2.12[d] displays the square envelop r 2 (f) and its 
estimate r 2 (t) in dB. Clearly, tracking is achieved between 0.4 ys and 0.5 ys, after the time 
0.5 ys, the estimator has a predicted value. 

Figure 4.2.13[a] shows the inphase and its least-square estimate. Obviously, the initial 
value of the inphase is 0.6, tracking is achieved between 0.8 ys and 0.9 ys; Figure 4.2.13[b] 
displays the quadrature and its least-square estimate. Clearly, the initial value of the quadra 1 
ture is 0.6, tracking is achieved during the time 0.8 ys and 0.9 ys; Figure 4.2.13[c] shows 
the mean square errors of the inphase and quadrature; Figure 4.2.13[d] displays the square 
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envelop r 2 (t) and its estimate r 2 (t ) in dB. Clearly, tracking is achieved between 0.8 ps and 
0.9 fj,s, after the time 0.9 fis, the estimator has a predicted value. 

Path #3, path #4, path #5 and path #6 illustrate the same properties as path #1 and 
path #2. Tracking is achieved during the different time interval for each path because there 
is a different time delay signal over it. 
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Figure 4.2.10: Zoom in Figure 4.2.9[d]: Path #1 - Path #3 



[a] Path #4, Time [u.s] 



Figure 4.2.11: Zoom in Figure 4.2.9[d]: Path #4 - Path #6 
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Figure 4.2.12: Multipath Component #1: Ricean Channel 
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Figure 4.2.13: Multipath Component #2: Rayleigh Channel 
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Figure 4.2.14: Multipath Component # 3: Rayleigh Channel 
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Figure 4.2.15: Multipath Component #4: Rayleigh Channel 
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Figure 4.2.16: Multipath Component #5: Rayleigh Channel 
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Figure 4.2.17: Multipath Component #6: Rayleigh Channel 
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Chapter 5 

Conclusion and Future Work 


5.1 Conclusion 

This thesis describes state space fading channel modeling and estimation using the 
stochastic differential equation. Models and estimates are obtained for inphase, quadrature 
and square envelope. 

1. Flat Slow Fading Channel, e.g., T s « T c . The state space estimates of I(t), Q(t), 
r 2 (t) show very good tracking at received signal-to-noise ratio as low and as -3 dB. 

2. Flat Fast Fading Channel, e.g., T s » T c . The state space estimates of I(t), Q(t), 
r 2 (t) also show that tracking is possible, when the received signal-to-noise ratio is 10 
dB. 

3. Frequency-Selective Slow Fading Channel. Similar results are also obtained for 
frequency-selective slow fading channels, which the state space estimates of the inphase, 
quadrature, square envelope of each path show very good tracking with respect to mean 
square errors, when the received signal-to-noise ratio is 20 dB. 

5.2 Future Work 

The various investigations pursued in this thesis, provide encouraging results, opening 
the possibility of continuing work in various directions. Some directions for future research 
are the following, 

1. Systems Identification. In this paper, all investigations have been done by assuming 
that the parameters of the state space channel models {A, B,C, D} are known. One 


67 



68 


CHAPTER 5. CONCLUSION AND FUTURE WORK 


can also estimate the channel states by first estimating the parameters {A, B,C, D} 
from measurement data. An algorithm for computing these matrices is given in [7] and 
[9]. 

2. High Order Approach. One can use high order rational transfer functions to ap¬ 
proximate the Doppler power spectral density function, and then employ order reduc¬ 
tion techniques to find lower order models which are closer to the true channel power 
spectral density. 

3. Robustness of the Models. The dynamic channel model can be described as a 
nominal one which its parameters are uncertain but having an range. Fixing the 
velocity v and the carrier frequency f c , then v/X is a constant. Thus, So(f) is a 
function of 0 m . Given a nominal transfer function H n (s) by choosing a 0 m . For 
appropriate selecting a weight function W 2 (s) and a variable stable transfer function 
A (s) which satisfies ||A(s)|joo < 1, one can get the perturbed transfer function which 
is defined by 


H p {s) = [l + A (s)W 2 (s) H n (s), s = ju. 


Based on the above perturbed transfer function, one can generate the robust state space 
channel model and then simulate the multiplicative perturbed inphase, quadrature and 
square envelope by using the robust Kalman estimation techniques. 

4. Receiver Design. Another interesting possibility for future work is to design receivers 
using optimal decision rules which do not assume knowledge of the channels. For the 
Bayes decision pay-off rule, the optimal decision rule can be computed explicitly by 
Mrosky to Kalman Filter estimate presented in this thesis. 

5. Optimal Coding Decoding. One of the most challenging problem is to design 
encoder and decoder which do not assume knowledge of the channel state information. 



Appendix A 


Expressions and Derivations 

A.l Notes about Aulin [1] 

A. 1.1 Part 1: Expressions 

Let 

E(t) = £ £„(<), 

71 — 1 

where E n (t) = Re{C n e^ n } at original coordinate, {C n , </>„, a n , P n }n=i are random variables, 
a n ~ u(0, 2n) and <f) n ~ u(0, 2i r). Here u(0, 2tx) is denoted the uniform distribution through¬ 
out 0 to 27T. 

(i) At any receiving point (x 0 ,yo,Zo)- Suppose that the transmitted signal is s(t) — 
Re{e^ ct }, then the received signal is 

E n (t) = Re{C n e j0J ^- Tn}+J ^}, 

where Re denotes the real part of the expression in its following brace, 

^c^~n 27T J'c^n 27 xf c 277”, 

c A 

A is the wave length associated with the carrier frequency, c is the velocity of light and d is 
the distance between the transmitter and receiver. That is, 

•<v 27T 

^ 7 '— Cn COST - d $ri) ■ 

A 

Let the vector (x 0 > yo,z 0 ) project to the n th wave path, then d = xo cos a n cos f3 n + 
y 0 sin a n cos f3 n + z 0 sin (3 n . Thus, 


69 



70 


APPENDIX A. EXPRESSIONS AND DERIVATIONS 


27T 

E n (t) = C n c,os[uj c t - — (x 0 cos cos P n + y 0 sin q„ cos (3 n + z 0 sin/? n ) + 0 n ], 

(ii) Let the point (x 0 ,yo,z 0 ) move in the received field. Then, (x 0 ,y 0 , z 0 ) —)■ (x 0 + 
vt cos 7 , y 0 + vt sin 7 , z 0 ), that is, 

E n (t) = cos [w c £ - ~ (^x 0 + vt cos 7 ) cos a n cos /3 n + (yo + vt sin 7 ) sin a n cos (3 n 
+ Zq sin fin'j + </>„] 

r 27T 

= C n cos u c t -— (vt cos 7 cos a n cos j3 n + vt sin 7 sin a n cos j3 n 

*■ A 

+ x 0 cos a n cos P n + y 0 sin a n cos P n + z 0 sin p n ) + (j) n 1 

n r + 2nvt 1 \ a 27rz 0 • a 

— On COS ]U c t -— COS(7 - OLn) COS /?„-— Sin p n 


2n Jx% + y% , 

--COS (a n + 0 COS Pn + (j)n\ 


— C n COS (u c t + U n t + On) 

— C n cos (u n t + 0 n ) cos u) c t - C n sin(o; n £ + 0 n ) sin u n t , 


where 


VJn ' 


2™ 

"X 

2ttz 0 

~T~ 


cos (7 - a n ) COS Pn 


2kJxI + yl 

sin p n - ■— y —, cos(o; n + 0 cos P n + p n , n e N. 


Here £ = tan 1 (yo/xo). Let 


N 


T c (t) = Yj Cn COS {u n t + 0 n ), 


n=l 

N 


Then we have 


T s (t) = Y c n sin(a ) n t + 0 n ). 

71=1 


E(t) = T c (t) cosu c t — T s (t) sintc c t. 


A. 1.2 Part 2: Mean and Variance 

(i) Take the expected value to T c (i). That is, 


N 


E[T c (t)] - E[Y Cn COS (u n t + 0 n )} 


n— 1 
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N 

= Y E[C n ]E[cos(u n t + 0„)], 

n=l 

since 

E[cos(u n t + 0 n )] = E[E[cos(uj n t + ^ n )U„=J] 

= E[~- / cos (ut + 0 n )dO n ] 

Ztt J o 

= o, 

then E[T c (t )] = 0. In the same way, we get E[T s (t)} = 0. Thus, 

E[E(t)\ = E[T c (t) cos u c t — T s (t) sin u c t] 

= E[T c (t)} cos u> c t - J5[T s (t)] sin u c t 
= 0 . 

(ii) Take the autocorrelation to E(t). That is, 

E[E(t)E(t + r)] — E[(T c (t) cosu c t - T s (t)sinu c t) 

x (T c (t + r) cos u) c (t + r) — T s (t + r) sino; c (t + r))] 
= E[T c (t)T c (t + t)] cos u c t cos u c (t + t) 

- E[T c (t)T a (t + t)] coso > c t sin a ; c (t + r) 

- E[T s (t)T c (t + t)] cos u c {t + r) sinu; c £ 

+ E[T a (t)T s (t + r)] s'uiUct sin u c (t + r), 

since 

N 

E[T c (t)T c (t + t)] = Y E\C n C m cos(u m t + 0 m ) cos(a > n (t + r) + Q n )\ 

n,ra=l 

l N 

= 2 E E[C n C m cos (co m t + oj n (t + r) + 0 m + 6 n )] 

n,m— 1 
1 N 

— A — "y ) E\C n C m COs{uJ m t Unit + t) + 0 m — @n)], 
^ n,m=1 

where 

{0 m + 0 n ) mod 27T ~ u( 0,27r), 

(0 m — 0 n ) mod 2-7T ~ u( 0, 27t). 
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Thus, we have 


E[C n C m cos(u J m t + Unit + r) + e m + 9 n )] = E[C n C m E{cos(u m t 

T LO n {t + t) + 6 m + 9n)\ 0 Jn,w m )] 
= 0 , 

E[C n C m cos (u m t - u n (t + r) + 6 m - 9 n )] = \ ^= l ^ C0SUnT ^ m = n . 

1 U, 771 z f = - 71/ 


Then, 


E[T c (t)T c (t + r)] = - Yj E [C n C m COS (U) m t - U n (t + T)+6 m - 9 n )] 

n,m—1 

= l Y E l C n cos u n r] 

Z n= 1 

= l Y E [Cl]E\coscj n T] 

Z n =1 

= RvY E i COSUJ nr}, 


where E[C^\ = Eq/N. Using the same method, we get 

E[T a (t)T,(t + t)] = ^Y E l cosu ^ 

n=1 

B[T e (t)T,(t + r)] = JjL X:S[sino,„r] 
£[T c (f + r)T,(t)l = -^E^Nn^r], 

71=1 


Thus, we have 


En " 


E[E(t)E(t + r)] = —— Y E[cosu n T] cos u c t cos u> c (t + r) 


n=l 

E ^ 

+ —— Y £-[cos uvr] sine J c t sin u) c (t + r) 

71=1 

E ^ 

+ — ~ Y E [ S i n UJ n T ] S ' n Wet COS U c (t + t) 
E n 

— — Y -^[sinavr] cosu; c isinu; c (t + r) 

n =i 
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where 


E n E ^ 

= ^ Yj E t cos u ”T] cos u) c r - ^ E [sin uj n r ] sin u c r 

n =1 n—1 

= E[T c (t)T c (t + t)] cos uj c t - E[T c (t)T s (t + r)] sina; c r 
= a(r) cos u c T — c(r) sin u c r , 


^ ^£[cOSUVr], 
n=l 

c(t) = ^^^[sinuv/-]. 


Let a ) n = tv, then we have 


c(t) = ^ ^E[sino; n T] 

71=1 

= ^^[sinwr] 

^0 nr • 1 

= —- E[smuT\ 
z 

— f sin(o;r )p u] du> 

E 0 f r 2n . , 2 -kvt 

2 J Ja—O A 

Eq r 1 . .27TUT 

— — —sin—-— 

2 J Je =o 27T A 

= 0 . 


cos (7 - a) cos (3)p a (a)pp((3)dad(3 


[ [ A sin(sing cos/ 3)pg(/3)d9df3 
J Je =o 27r A 


= ^E£[ cos Wn r ] 

= f^ NE i COSU,r \ 

= ^-E[cosu>t\ 

— ^ J cos (urfpudu 

Eq ft 2 * , 2lXVT 


cos (7 — a) cos P)p a (a)pp((3)dad(3 


n fZTT / 7 T 77 T* 

/ / cos(-—-cos( 7 - a)cos(3)p a {a)pp((3)dad(3 

Z J J Q!=0 A 

T / L- o i cos (^Y- sin (f - (7 - <*)) cos (3)pp(p)dadp. 


Let 6 = 7 r /2 — (7 — a), then 


a( r ) = ^ f [ 2 Acos(^y^- sindcosP)pp(P)d,QdP 

2 J Je=?r—i 2n A 
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Eq 

Y 

Eq 

Y 


J -^-pp(f3) J cos^^Y- sin9 cosP)d9d(3 

j cos (3)p p ((3)dp, 


where Jo(-) is a Bessel function of the first kind of zero order. That is, 

Jo(z) = A- f e jzsine d9. 

Ztt J —7r 

A. 1.3 Part 3: Rayleigh Fading 

(i) Suppose pp(/3) — 5(f3), in this case, denote a(r) by a 0 (r), then 

a 0 (r) = a(r) 

= Y J cos(3)8(f3)dp 

E 0 2twt 

- T m ~>’ 

Take the Fourier transform to a 0 (r), we get 


Mf) = E[a 0 (r)} 

1 E 0 r°° r e j2TTVTsine/X e -j2iTfr 

J—oo J—n 


2 tt 2 

F 0 

r 

47T J 

e= 

Fo 

r 

47t J 

e= 


Let / = | sin#, then 9 = arcsin | IE | = s i n 0 < p Thus, we get 
d9 = 

Finally, we obtain 


, 1 - df , | / |< L »€ 

NEW” A 


7T 7T 
2 ’ 2 


A 0 (/) = ir- f S(f-Y'm9)d6 

47T J0=— 7t A 

= ^ / X £(/ — y sin#) 1 == -d/ 
4tT Jf=-l A yh _ (XA)2 V 


Ec\ A 


(ii) Suppose 


4 "VhW 

1 0, otherwise 
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f N Eq 1 fPm ZTTVT 

a CO = V 9g - cos /3) cos (3d(d. 

2 2 sin Pm Jb — 8m A 


2ttvt 


2 2 sin PmJp^-Pm 
Take the Fourier transform to a(r), we get 1 


A(f) = T[a{r)\ 

= — — — [ 0m [" [°° cos pe j2 ™ TCOSl3sine/x e- j27TfT dTd9dp 

8 tt sin /3 m Jp=-p m Je=-n Jt =-oo 

En fP™ /' 7r r> 

= -— / / 8( f ——cos/3 sin 9) cos fid9d/3 

8ir sm p m J/ 3 =-/? m Jfl=-7r A 

En [Pm rit/1 y 

— -—;—— / / <5(/ — — cos/3 sin 0) cos fidddft. 

2it sm /3 m Js=o Je=-n /2 A 


27t sin /3 m Jj3= o 


x = cos (3 

y = v^ sin0 


j sin/3 = 

0 T 

- —x^ I 

-\- =4> S 

f M < 

4 

„ 7T 7T1 

=> 0 e - - 

[ cos 0 = 


_ X 2 ! 

*u v 1 

VI 

S& 

_ / 

4 

2 2-1 


Since Jacobian determinate is 


<*(s,y) = -]/i sin P 

d(P, 9) 0 


^ cos 9 


’i-Ni-y'- 


Then, 




27T sin f3 m Jx^yf^cospm ''V=-yf% — X “ y 2 

< 0 , 


I/I > h 


■So 1 A 
2 2 sin p m v ’ 


Eg 1 X I n 
2 27r sin P m v \ 2 


f cos Pm 

< I/I <!; 


- arcsin 2cos , |/| < Jcos /3 n 


1 As the integration is symmetry, cos a is an even function. Thus, 

n n n/ 2 


n n 

cos ada = 2 / cos ada 

-7 Jo 


-jC 


cos ada. 
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(iii) The envelope R(t ) and the phase <f>(t) are defined as follows 

f T c (t ) = i?(f) cos<f>(t) 

\ T a (t) = R(t) sin §(t) ’ 

then R 2 (t) = T 2 (t) + T 2 (f). Since E[T c (t)} = E[T s (t )] = 0, E[T e (t)T c (t + r)] = a(r), 
E[T s (t)T s (t + r)] = a(r), thus the second moment (variance) a(0) = jB[T c 2 (t)] = -E[T s 2 (t)]. 

Let y = T 2 (t) + T 2 (t), then y is chi-square distribution with two degrees of freedom. The 
pdf is 

'w-slo r ,/M0> - yi0 - 

Then 

Let R(t ) = + T 2 (t) = y/y, then y — R 2 (t), dy = 2 R(t)dR(t). Integration the pdf 

with respect to y gives 

r p ^ - r wf* /i ‘ { ° )dR =l 

Thus, the pdf of Rayleigh distribution is 


PR(r) = ^) e * rV2 “ (0> ’ r - °- 
And the mean of the envelope is 

f OO 

£[£(*)] = l rP R (r)dr 

= r°4^* 

7o a(0) 

/•OO 

= f -re~ r /Mo) d(-r 2 /2a(0)) 

Jo 

+ r e - r2 M° Mr 

Jo 


= —re 


—r 2 /2a(0) 


/ e - r2 / 2a (°)^ 
Jo 

\/2tto(0)^ 

IfaaiO) 


o 

'r 

oo ] 


y27ra(0) 


2 /2a(0) dr 


7ra(0) 


2 



A. 2. A UTOCORRELATION FUNCTION AND PSD 


77 


A. 1.4 Part 4: Ricean Fading 

For the specular component with the parameters {a Q ,p 0 , <f> 0 , \REd}, we have the envelope 
R r (t), which is Ricean distribution 2 , 


PRriT) 


o> /o(r gp, r > 0 


where Iq(-) is the zeroth-order modified Bessel function of the first kind. 


A.2 Autocorrelation Function and PSD 

(i) Consider a wide-sense stationary input signal x(t ) which passing through a linear 
time-invariant system having impulse response h(t). The output signal y(t) is given by 

/ OO 

h(r)x(t — r)dr 

-OO 

in time domain, which has the format 

Y(f) = H(f)XU) 

in frequency domain. Here Y (/), X (/), //(/) are the Fourier transforms of the output signal, 
input signal, impulse response. The autocorrelation is given by 

Ryyir) = E[y(t + T)y*(t)\ 

/ OO p OO 

h(p)x{t + T- (3)d[3 / - 7 )d 7 

-oo J —OO 

= / / h*('y)h(p)E[x*(t-^f)x(t + r-p)]djdp 

/ oo poo 

/ h*(l)h(P)R xx (r + 7 - p)dnfd/3 

-OO J — oo 

= J ^ (■Rxx('r + 7)=t=fi(r + 7))fi*(7)d7. 

Let 7 — —£, then 

= J_^(Rxx{r + 'y)*h(r + j)^h*('y)d'Y 

= J (R xx {t -t)*h(r -t)}h*(-t)dt 

= Rxx{t) * h(r) * JT(—t). 

2 there are two ways to find the pdf, one is on page 43 of [17], the other is like the same way of finding 
the Rayleigh. 
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(ii) Take the Fourier transform to Ry y (r), we obtain 
s,(f) = r R„(r)e-^dr 

J —CO 

/ oo poo poo 

/ / h*(7)/i( / 5)72 ;EX (r + j-P)e~ j2nfT d'ydpdr. 

-oo J—oo J—oo 

Let u — t + 7 — /?, then 

/ OO /-OO r OO 

/ / h'(j)h(l3)R.*(r+')-0)e- 12 ’ ,T d'ld0dT 

-oo J —oo J —oo 

= /°° /°° /°° h\7)h{P)R xx (u)e-^ u+ ^d^dpdu 

J —oo J—oo J —oo 

/ OO /'OO /'OO 

/ l *(_ 7 ) e -i27 r / ( - 7 ) d (_ 7 ) / h{(3)e~ j27Tf/3 d(3 / R xx {u)e- j2nfu du 


A.3 Approximation DPSD by // s 

A. 3.1 Approximation 71 ( 5 ) 

Let the second-order rational transfer function be 

k 

s 1 + 2Cw n s + o;2 ’ 

where s = jfu;. Then 

if CM = -——. 

- w 2 + 2(ju n u 

Thus 

I I = \jH(ju)H(-ju) 

_ k _ 

yj ( 0/2 - a / 2 ) 2 + 4 C 2 w 2 o ;2 

Define 

fmax ,/ci COs(/3 m ), 


(A.3.1) 

(A.3.2) 


« = 
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that is, 
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because the property of the damping ratio ( -» 0 as Soifmax) —> Finally, we have the 
parameters 




Sp( 0) 

SD(fmax) ’ 


0J r 




and k — uI^Sd(0 ) 


for the transfer function H(s). 


A.3.2 Multiplicative Perturbation of H(s) 

Replace Soif) by A(f) in the previous section, we got the parameters {(, k, u n } of H(s) 
for the normalized DPSD, which only depends on /3 m . 

Let the nominal transfer function be 

H n(s) = ,, ■ ■ 

S CLfiS “1” Ufi 

and the perturbed transfer function be 

H r (s) = ^ 

s 2 o,pS -t~ bp 

then we have, by the definition, the weight function 


W 2 (s) = 


H P (s) 

H n {s) 

(fe ~ 1 ) s2 + (fe Q » ~ a p) s + ~ M 

S 2 ~b CLpS bp 


A.4 From H(s) to State Space Representation 


Consider the Laplace transform of the impulse response 




Y(s) 


k 


U(s ) S 2 + 2((jJ n S + u> 2 ' 

(i) The Controllable Canonical Form. Based on the equation above, we define 


*00 

Y(s) 


s 2 + 2(uj n s + 
X(s) 


— feed back, 

— feed forward, 


then we have, by taking the inverse Laplace, the following equations 
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£ _1 {(s 2 + 2C u n s + oJ 2 n )X(s)} = ^{[/(s)}, 

£- 1 {y( s )} = £- 1 {x( s )}. 

From (A.4.3), we get 


x + 2(u n x + u? n x = ku. 


Let zi = x, z 2 = x, that is, ii = z 2 , z 2 = x, then we have 


z\ = z 2 

z 2 = ku - 2(u n z 2 - uP n Z\ 


that is, 


ii 


0 

1 

-Zi 


'o' 

-^2 _ 



2C ) D n 

/2_ 

+ 

k 


from (AAA), we obtain 


(AA.3) 

(AAA) 


y = x - zi 


that is, 


y = [ 1 


0 ] 


zi 

Z2 


(ii) The Observable Canonical Form. From the equation 


Y(s) _ k 

U (s) s 2 + 2 (u n s + u; 2 ’ 

we get 

Y(s) = -^Y(s) - Ay (s) + 4u(s) 

= )( - 2Cu„y(.s) 4 U- t,£Y(s) + kU(s)>). 


ZiM = )(-ulY(s) + kU(s)) 
Z 2 (s) = l(-2(u,„Y(s) + Zt(s)) ’ 


(AA.5) 


then 


Y(s) = Z 2 (s). 


(AA.6) 
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By substituting equation (A.4.6) into equation (A.4.5) and multiplying both sides of the 
equations by s, we obtain 


f sZ 1 (s) = -ulZ 2 (s) + kU(s) 

\ sZ 2 (s ) = — 2(u n Z 2 (s) + Zi(s) 

Take the inverse Laplace transforms of the preceding two equations, we get 


that is, 


and 


zi(t) — -u1z 2 {t) + ku(t ) 
Z 2 (t) = — 2^co n z 2 (t) + Zi(t) ’ 


Z\ 

A 

= 

0 

1 

“Wn 
—2 (u n 

Z\ 

z 2 

+ 

1-1 

^ o 


y = [o i' 


Zl 

Z2 


A.5 The Kalman Estimator 

A.5.1 The Noise Covariance Matrix 

Let the additive bandpass noise is v(t) = v/(t) cos(cu c t) — vq(I) sin(u; c t), {vi(t)} t > 0 and 
{vQ(t)}t>o are f wo hd with density N(0; a„). Then the covariance matrix is 

R = E[v(t)v tr (t)} 

E[vj{t)} Elv^VQit)} 

E[vi(t)v Q (t)) E[vq(t)] 

‘ E[vj(t )] 0 1 

- 0 E[v* Q (t)]\ 

'*1 o ' 

0 a 2 v ’ 


and so does the Qi. 

" We know the state space model 


x — 


y 


Ax + Bu 
Cx + v 
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with error covariance equation and Kalman gain. We want to find the error covariance 
equation and Kalman gain of the state space model as follows 

x = Ax + Bu 
y — Cx + Dv 

Let the state space model be 

x — Ax + Bu 
y = Cx + v' 


where v' = Dv, then we get 


R! = E[v'v ,tT } 

= E[Dvv tr D tT ] 
= DE[vv tr ]D tr 

= “L 

L 0 °v\ 

= DRD tr . 


A.5.2 Ito Method 

Consider a time interval [0, T] with 0 = t 0 < 4 < • • • < 4 < • • • < £jv = T, which has 
step size 5 k = t k + 1 — t k . For equation 

dx = Axdt + Bdw, 

we have 

x(t k+ i) - x(t k ) = Ax(t k )(t k+1 - 4) + B(w(t k+ i) - w(t k )), 
where x € Let Aw k = te(4+i) — w(t k ), then we get, by denoting k = t k , 

x k+ i = Ax k 6 k + x k + BAw k 
= (A5 k + I)x k + BAw k , 

where E[Aw k ] = 0, E[{Aw k ) 2 ] = 6 k , I e R n x Sfi” is an unit matrix. Let S k — T/N — 5, 
which means that the step size has an equal length. Denote Aw — Aw k , then E[Aw\ — 0 
and E[(Aw) 2 ] = 5. That is, Aw ~ N(0]5). Thus, 

x k+1 — Ax k 5 + x k + BAw 
= ( A5 + I)x k + BAw 
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A.5.3 Computation of the Conditional Distributions 


We know that if the vector of random variables ( X , Y, Z ) is joint Gaussian with zero 
means and covariance matrix 


Mxyz = 


var(X) cov(X,Y ) cov(X, Z) 
cov(Y,X) var(Y) cov(Y,Z) , 
cov(Z,X) cov(Z,Y) var(Z) 


then the marginal pdf for X and Y is also Gaussian and has the same set of means, variances 
and covariances 3 . That is, the pair (X,Y) has zero means and covariance matrix 


M xy = 


var(X ) 
cov(Y,X) 


cov(X, Y) 
var(Y) 


(i) Conditional Envelope Distribution. From the least-square estimation of the 
conditional stochastic process, we have 4 


I(t) = E[I(t)\Yt] and Q{t) = E[Q(t)\Y*\, 

where Yq = (y(s); 0 < s < t}. Then the conditional joint Gaussian pdf of the two stochastic 
processes, the inphase and quadrature, is 




^ -b(x-x) tr M 1 (x-x) 


where 


x = 


■ m ■ 


' At) ' 

Q(t) 

5 X 

Q(t) 


Next, we will find the pdf of the envelope given by Yq. Let Z = I 2 + Q 2 , then we have, 
by the definition, 


Prob(Z < z\Yq) - Prob(I 2 + Q 2 < z\Y£) 

= f f(I,Q\Y*)dIdQ 

JI 2 +Q 2 <z 


-l 


I 2 +Q 2 <z 


e -±(x-x) tr M 1 (.*-*) dldQ. 


Let the covariance matrix be 


3 see [15] 
4 see [13] 
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M = 


a b 
b c 


then the inverse matrix of M is 


M - 1 = 


Thus, we have 


1 


ac 


-V 


c —b 
—b a 


--(x — x) tr M 1 (x — x) = 


g rn 
m n 


1 

~2 

\LJ L. -I / L. -J\ 

-§(/ - Tf - m(l - /)(<? - Q) - ~(Q - Qf 

-f P-^Q 2 - mlQ - 9 -P - \Q 2 - mIQ 
+ (gl + mQ)I + (nQ + mf)Q. 


/ 

/ 


T 

\ tr 

5 

m 




l 

<9 


LQ. 

) 

m 

n 

l 




Let 


(I = v /pcosO 
10 = 


Q = y^psin^ ’ 
then we get 0 < p < z, 0 < 0 < 2it, \ J\ = and 

r 2n 1 


n Z7T I 
2 


with 


~(x — x) tr M 1 (x — x) = —^T 1 — r ~Q 2 — mIQ 

9 „ , ( 9 -n)p -2 


~ 2 P + 


sin 2 0 — mp sin 0 cos 6 


+ (gl + mQ)^fp cos 9 + (nQ + ml)y/p sin 6. 


Then 


p 27 T 


n Z7T I 

-f(p,e\Yt)d P de 

= fl- 

■'o 2 2 tt v '|M| - 70 

x g(nQ+rn/) ^/psinQ & p sin 2 9 & -mp sin 9 cos 6> ^g 


(gl+mQ) y fp cos 6 


<o> 
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Since e x = ~r> we have 5 

f g (gI+mQ) y/p cos 0 ( nQ+ml) y/ps'md^ ps'm 2 6 —mp sin 9 cos 6 

Jo 


°° (—V\ k 2~i ^ 

= 4 £ hUUuA ( gI + m Q) h ( n Q + mI ) l (9 - n) 3 m k 

h,i,j,k =0 

h+i+2 i+ 2k [* „h+k 


Xp 2 

oo 


j : 


cos h+k $ shT +2j+fc Odd 


= 2 Y, Lll (9l + mQ) h {nQ + ml)\g - n) j m k 


xp 


h,i,j,k =0 

fa+i+2 , +2fc r(^p)r(^ +2 y -hi) 


p^ ^ , ~t~^~^2j-{- , 2fc-|-2 ^ 

Thus, we have 

Prob{Z < z\Y<) = [~ - l—e-P-ffi-nW-b 

J o 2 tx^\M 


00 /_i ^ ^ ^ 

x £ 71 -i..., {gl + mQ) h {nQ + ml)\g - n) 3 m k 

h,ij,k =0 n - l -3- K - 

^±±±2£±2fc -L 2 ( 2 / 


xp 


p^ fo-H~h2.7'4~2AH~2 ^ 

Then the conditional pdf of the channel square envelope is 


-dp. 


fa{z,t\Y*) = 


1 _ &-mIQ-%i 


2 TTy/\M 

~ (—l) fc 2“ J 

X h ,^o~MW 


(gl + mQ) h (nQ + ml) l (g — n) j m h 


x 


r(fe±|±i)r(i±2T+T+i) „ +i+2j+2fc 


p ^ +2fc+2 ^ 

Let r = s/z = a // 2 + Q 2 , then z = r 2 , dz = 2rdr. Thus, we obtain the conditional pdf of 
the channel envelope 

f(rtV «) 1 .-«P-*0»-mES-*,* y- (-1)‘2-U(4±|±l)r(tta±4±l) 

/a( ’ 1 o) - AKWIftl r(4±i±^) 


x((/J + mQ) h (nQ + mT) l (g — n) j m k r h+l+2j+2k+1 . 

(ii) Conditional phase distribution. Let 6 = tan -1 (^), then 
Pro&(0 < 0|T O *) = Pro 6 (tan _1 (Q\ < <f>\Y$) 


5 see [ 16 ] 
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Let 


= L^)<j ( ™ IdQ 

= [ . -- \=e-^ x -^ trM ^ x -^dIdQ. 

./tan -1 (y )<</> 27ry |M| 

I = p cos 9 


{ Q = psin0 ’ 
then 0 < </>, 0 < p < oo, | J| = p. Thus, we have 


Prob(9 < (/>\Yq) = [°° [ 

Jo Jo 

= n 


oo r<f> \ 


where 


2?r^/|M 


e -i(x-S)‘r M \x-Z) pdedp 


[* 1 c -jfi-gQ 2 -mJQ-gp 2 

2^ 

x e (p/+mQ)pcos0 e (nQ+m/)psin 0 g sin 2 6 e ~mp 2 sin 0 cos ® pd Q dp 

= f°° [* - 

■/o 'o 27ry|M 


J? = 


19 9~n„. r 2 

2 2 


sin 2 0 + m sin # cos 


Thus 6 , 


5 = 


(p/ + mQ) cos # + (nQ + ml) sin 0 
2P ‘ 


Pro&(0 < </>|Y 0 *) = [* - 1 &- mI Q e s2 d9 

JO O-rrJlMl 


X 


b 

l 

2lXy 

Jw\ 

roc 

' -1 

k 

2 P 2 ° 

b 

1 

- P. 


S 

P* 


27r\/iMi 

Thus, the conditional phase distribution is 


y/2 R' 


MY$) = 


A _ 1 + Sy/ne ^_ap_ 2 Q 2 _ m/ Q 


4 7 rP 2 y / |M 


-e 2 2 ' 


6 /o°° e £*dx = ^a, because = 1 


+°o 1 
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A.6 Ricean Channel Model 

A.6.1 State Space Representation 

Let E[I R (t)\ = ro cos(uj()t + <fr 0 ), then the Ricean state space channel model can be 
represented as follows 


X R {t) = Xj(t) + r 0 cos(a j 0 t + fro) 


(A.6.7) 


where Xi(t ) is Rayleigh case, the parameters {r 0 ,o;o, fro} are given. Take the derivation to 
(A.6.7), that is, 


Xp(t) = X T (t) - r 0 uJo sin(w 0 t + fro) 


A I X I (t) + Biwi(t) - r 0 u > 0 sin(o;ot + fr o) 

- r 0 cos(a> 0 t + fro) q ) + BiWj(t) - r 0 u 0 si n(u> 0 t + fr 0 ) 


tqoJq sin (iv 0 t + (fro) 


= AjXf (t) + - AiTq cos(u Q t + <fro) 

= A.Xfit) + Bjwjtt) + 

Let E[Q R (t )] = ro sin (u 0 t + fr 0 ), then as the same way of finding X R (t), we get 


-r 0 u) 0 sm(u) Q t + (fro) 
r 0 ul cos(uj Q t + (fro) 


X R (t) = A Q Xg(t) + B Q w Q (t) + 


rR 


vqujo cos(ujot P fro) 
r 0 ul sin(o;ot + (fro) 


Thus, define A 
given by 


a r At 0 1 „ A 


1 0 A, 


Q 


,B = 


Bi 0 

0 B Q \ 


, the state space representation of {y(t)}t >o is 


X R (t) = AX R (t) P Los{t) + Bw(t) 
y{t) = G(t)X R (t) P D(t)v(t) 


where G(t) = s(t)C(t), 


-r 0 uJo sin (u Q t P (fro) ' 
r Q u)l cos(tu 0 t P (fro) 
r 0 uo cos (u 0 t P (fro) 
roufr sin(w 0 t + fro) . 


Los(t) — 
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A.6.2 The Filter for the Ricean Channel Model 


Since the Ricean and Rayleigh have the following relationship 


X R = X + 


r 0 cos(u 0 t + 0o) 
0 

r 0 sm(u 0 t + 0 0 ) 

0 


Take the derivation to the above equation, we get 


j? 


X = x + 


-r Q u 0 sin (u 0 t + 0 O ) 

0 

r 0 u;o cos (u Q t + 0 O ) 

0 


AX + K(y - GX) + 


AX + K(GX + Dv- GX) + 


-r 0 uJo sin(a j 0 t + 0 O ) 

0 

ro^o cos (u 0 t + 0o) 

0 

-r 0 u> 0 sm(u 0 t + 0o) 
0 

tqujq cos(tao t + 0o ) 
0 


r 0 cos (u Q t + 0 O ) 

0 

r 0 sin(u! 0 t + 0 O ) 


= a(x r 

L 

+ Dv - G(X r - 
= AX r + k(GX r + Dv-GX r ) 
- A 

= AX H + Los(f) + if(?/-GX K ). 


) + k{g(x r 


r 0 cos(u 0 t + 0o) 
0 

r 0 sin(u; 0 £ + 0 O ) 


0 J 


L o j 

' r 0 cos (u 0 t + 0o) ' 


—rou>o sin(u 0 t + 0o) 

0 

)) + 

0 

r 0 sm(u 0 t + 0o) 

r ou»o cos(a;ot + 0o) 

0 


0 


' r 0 cos (u 0 t + 0o) ‘ 


' -r 0 oJo sin (u 0 t + 0o) ' 

0 


0 

r 0 sin((u 0 t + 0 O ) 

+ 

rou>o cos(uj 0 t + 00 ) 

o 


0 
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A.7 Signal-to-Noise Ratio 

A.7.1 White Gaussian Noise X(t) ~ N(0]No) 

(i) The definition of the power spectral density function 7 is 

Sx(f) = HRx( r)} 

- 

Rxir)e-^ T dr, 

-OO 

where X(t) is a continuous time WSS random process with mean rrix and autocorrelation 
function R x (r), which defined as follows 

R x (r) = E[X(t)X(t + r)} 

= N 0 S(r). 

Especially, Rx( 0) = E[X‘ 2 (t)] — Nq. Thus, we have 

/ OO 

Rx(T)e-*' tr dT 

-OO 

/ OO 

N 0 5(T)e- j2nfT dr 

-OO 

= No, 

which means that if X(t) ~ N(0, a 2 ), then the power spectral density of X(t) is a 2 . Let the 
bandpass additive noise (substitutes the user’s interference, but it is the worst case because 
the user’s interference has band limit) is 

v(t) = vj(t ) cosu c t — vq{P) sine o c t, 

where {f/(i)}t>o, {vQ(t)} t >o are two iid white Gaussian noises with zero-mean and variance 
N 0 . Then, 

E[v 2 (t)] = E[(vi(t) cosu c t — VQ(t) sinu c t) 2 ] 

— E[vj(t)\ cos 2 a } c t + E[vQ(t)\ sin 2 u c t 

= E i v i{t)] 

= N 0 . 

(ii) Let ( SNR t )dB be the transmitted signal-to-noise ratio in dB, E t be the transmitted 
signal power. Then, 

7 see P404 of [15] 
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1 r T s 2 

Et / I ^transmitted signalil) | dt, 

-l s **'0 


where T s is the signaling interval; and 


{SNR,) dB = 101og 10 §-. 

iv 0 


Thus, 


iVo = 




(SiVR t ) dB 

10 10 


which means if we have white Gaussian noise n ~ JV(0,1), then 

E t 


E, 


1 0 
0 1 


A , n ( SNR t)dR n — v i ~ N(Q,N 0 ), R— (SNRt)dR 

\| 10 io 10 io 

(iii) Let ( SNR r )dB be the received signal-to-noise ratio in dB , E r be the noiseless received 
signal power. Then, 


1 f Ts 2 

E r / | S rece i ve d signal (0 | dt. 

1 c JO 


Thus, 


N 0 = 


E r 


(SNR r ) dB 
10 10 


which means the additive noises are 


vi 


E r 

ni 

V Q_ 


(SNRrUg 
10 10 

n Q 


where {n/,(£)} t >o ~ N( 0; 1), L — I,Q. And the covariance matrix is 

E r 


Rv = 


( SNR r ) dB 

10 10 


1 0 
0 1 


A.7.2 The Noise Power 


By the definition, the average noise power is 

rT, 


K = ~E[J^ \D(t)v(t)\ 2 dt] 

[ T ‘ D(t)E[v(t)v‘ r {t)\D tr (t)dt 
Jo 


l r Ts 
T s 
1 r T ■ 


i r J * 

“ T,L d V 

= No. 


E[vj(t)) 0 

0 E[vl(t)}\ 


D tr (t)dt 
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A.7.3 The Rayleigh Channel Received Signal Power 

By the definition, we have 

Er = Le[ P | G(t)X(t) 

J s ^ 0 


dt 


= If 

Is Jo 


E 


G(t)X(t) 


dt 


- f Q S s 2 (t)C(t)E[x(t)X tr (t)]C tr (t)dt 

rp 

J q 3 sW^EiXj.X*] cos 2 (u c t) + E[X Ql X% J sin 2 {u c t))dt 
[ S s\t)E[X h X%\dt 

T 

[ s 2 (t)7n(t)dt. 

Jo 


1 rT ‘ 
T. 

1 * 
T s 

1 fTi 
T s 


aJ Jy al+OIf {s(t)} t >o is a narrow band signal, and hence without loss of generality 

assumed to be s(t ) = 1, by performing the above integration, we get 

E — 2C 2 A — 2 (3a ~ Ps ~ Pi r -iCu n T s 


+ 


4 (u n T s “ +A 

~ C 2 + Wi + 2C 2 /?5 + C 2 /^3 + /^2 — C 2 /^2 + /?3 
4 (AVX 


+ - 2 ^ 5 - 2 flv /r3 7 cos (^A Z V T -) e ~ W - 

4u) n T s v v / 

+ sin (2 u , n JEP2 T \e-^- T - , 

4w n T s v v ' 

where 7 (0) = £pf(0)2sT <r (0)], 

n ( 722(0) +711(0) Cw„ k 2 \ _ 71l(0)C 2 ^n+ 2 ')'12(0)CWn+722(0) 

40,3^/HciJ’ /J 2 “ • ^(1-C 2 ) 


P 3 = 7ii(0), 






4u>3(l-C 2 ) 1 


fio — 


4C^3(1-C 2 )' 

k 2 


A.7.4 The Ricean Channel Received Signal Power 

Let s(t) = 1, then we have, by the definition, 

E? = Ie\ f \ G(t)X{t) + G(t)h,(t) 


dt 



A.7. SIGNAL-TO-NOISE RATIO 


93 


l r T > 


f 


T t 
l r T ‘ 


E 


| G(t)X{t) + G{t)h s (t ) | 2 


dt 


- J o S s\t)C(t)E[x(t)X*(t)]c*(t)dt + yI° 

TjT J o 1n(t)dt + J-J 0 C(t)E[h,(t)h‘;(tj}C ,r (t)dt 


±r 

T s Jo 
— E r -f- los , 


where 


h s (t) = 


r 0 cos (u 0 t + 4> 0 ) 
0 

To sin(ui 0 £ + (p 0 ) 
0 


1 T 

los = -jf ' C(t)E[H,(t)H?(t)]C‘ r (t)dt 

= ~ f (tq cos 2 (oi c t) cos 2 (u) 0 t + 0o) - 2ro sin(u; c t) cos(cj c t) 

J- s ** 0 

x sin(n> 0 £ + cj) 0 ) cos(<jj 0 ^ + </>o) + r% sin 2 (uj c t) sin 2 (uj Q t + <f>o))dt 
r 2 /-r s 

= -f / cos 2 (cu c t + u) 0 t + fa)dt 
j s ** 0 


rg r o s ^ n ( 2 (w c + a, 0 )T s + </>o) 
2 4(a> c + u>o)T s 


A. 7.5 The Multipath Received Signal Power 


E, 


l ~ E \ [ Tm \J2GiX i \ 2 dt 

m JO i=1 

= E\{Y,G t X t )C£G i X i ) 

m J0 L i=l j =1 

= E E\Y,G i X l xrGT\dt+ ^ G E[ £ G.X.Xfff 
m Jo L »=i Jm Jo 1,7=1 


dt 
1 f Tm 


dt. 


*,3= 


Since 




E. 


S{t - Ti)s{t - Tj) = | 1<*,J<6 


1 r rT m 6 

—E 


7 m i x>*< r * 

L ‘ /0 7=1 




Then 
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1 p'JT 6 -t rp 0 

= 7T / " E\'£G i X i X?G?]dt + - r f m E\t G i X i Xj r Gf]dt 

Im J 0 T m JO i,j=l 

= A f T "’£G,E[X l X!']G^it 

J- m ^0 

= [ ^s 2 (t — Ti) (7J1 + 
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